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20.  ABSTRACT  (Continued) 

^ explained  by  the  release  of  reasonable  amounts  of  neutral  vapor.  The 
subsequent  deposition  of  ionization  in  the  expanding  nOutral  cloud  results  in 
an  ion  cloud  that  is  substantially  narrower  than  the  neutral  cloud  consistent 
with  a number  of  radar  measurements.  The  influence  of  neutral  winds  and 
settling  under  gravity  on  the  distribution  of  ionization  along  magnetic  field 
lines  and  the  descent  through  the  atmosphere  is  treated  analytically.  It  is  ^ 

shown  that  downward  winds  can  delay  the  decrease  in  the  peak  ion  concen-  » 

tration  and  that  altitude-dependent  winds  are  more  effective  in  limiting  the 
increase  in  length  of  the  ion  cloud.|_A  quantitative  model  of  the  Spruce  ion 
cloud  based  on  these  new  modelingq>^cedures  is  developed.  The  best 
estimate  of  the  total  ion  inventory  is^±  1 x 1024.  jhe  transverse  scale 
of  the  ion  cloud  is  2.66  km  and  the  tot^  ion  content  on  the  field  line  passing 
through  the  center  of  the  ion  cloud  is  2^  x 10  m“2.  The  peak  electron 

density  of  2.45  x 10^2  ^-3  is  reached  aAapproximately  30  seconds  aftei' 
release  and  downward  neutral  winds  and/or  decreasing  atmospheric  density 
scale  heights  serve  to  keep  the  peak  ion  cmcentration  of  the  ion  cloud  above 
7 X 1012  m" 3 for  more  than  30  minutes.  \ 
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SUMMARY 

We  have  developed  a model  that  quantitatively  describes  properties 
of  ion  clouds  produced  from  barium  releases  in  the  lower  ionosphere.  First, 
a quantitative  procedure  is  described  for  the  formation  of  the  ion  cloud  from 
the  expanding  cloud  of  neutral  barium  vapor.  Second,  the  subsequent  expan- 
sion of  the  ion  cloud  along  magnetic  field  lines  and  its  descent  through  the 
atmosphere  under  the  influence  of  gravity  and  neutral  winds  is  treated  analyt- 
ically. Third,  the  results  of  these  procedures  are  applied  towards  determin- 
ing quantitative  characteristics  of  the  Spruce  ion  cloud. 

A description  of  the  early -time  expansion  of  the  neutral  barium -atom 
vapor  cloud  is  described  by  a snowplow  model  with  diffusion  added.  The  addi- 
tion of  diffusion  aids  in  making  the  transition  from  an  early  rapidly -expanding 
dense  cloud  to  a late -time  cloud  that  is  expanding  solely  by  diffusing  through 
the  ambient  atmosphere.  It  is  shown  that  the  effect  of  diffusion  is  to  rapidly 
change  an  arbitrary  initial  density  distribution  into  a Gaussian  distribution 
which  then  begins  to  increase  the  radius  of  the  neutral  cloud. 

A detailed  examination  of  the  results  of  the  theoretical  model  in 
comparison  with  a film  record  of  the  Spruce  neutral  cloud  at  20  seconds  after 
release  is  made.  This  comparison  provides  a plausible  explanation  for  the 
observed  fact  that  photogrsqihic  images  of  early -time  neutral  clouds  appear  to 
be  far  larger  than  could  readily  be  explained  by  the  release  of  reasonable 
amounts  of  barium  vapor.  Previous  modeling  attempts  have  required  the 
assumption  that  more  than  twice  the  amount  of  vapor  was  released  than  is 
assumed  in  this  report. 
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The  subsequent  deposition  of  the  ion  cloud  from  the  expanding  neutral 
vapor  cloud  is  treated  quantitatively.  The  results  show  the  ion  cloud  is  sub- 
stantially narrower  than  the  neutral  cloud.  The  result  of  this  model  deposition 
is  to  produce  an  ion  cloud  with  transverse  widths  which  are  in  substantial 
agreement  with  radar  measurements  of  the  transverse  scale  size  of  the  early - 
time  ion  cloud.  We  find  that  the  high -density  central  portion  of  the  ion  cloud 
is  more  peaked  than  a Gaussian  distribution,  but  that  the  outer  regions  contain 
far  more  ionization  than  a Gaussian  contains. 

The  subsequent  development  of  the  ion  cloud  along  magnetic  field 
lines  and  its  descent  through  the  ambient  atmosphere  is  treated  analytically . 

The  treatment  includes  a magnetic  field  of  arbitrary  dip  angle,  the  presence 
of  neutral  winds,  and  examines  the  responses  of  both  low-  and  high-conductivity 
ion  clouds.  It  is  found  that  the  time  to  assume  an  asymptotic  shape  is  sub- 
stantially longer  than  had  been  estimated  previously.  Downward  neutral  winds 
can' substantially  delay  the  decrease  in  the  peak  ion  density  and  the  spread  of 
the  ion  cloud.  Altitude -dependent  neutral  winds  appear  to  have  a greater  effect 
in  limiting  the  increase  in  length  of  the  ion  cloud.  Moderate  upward  winds  do 
not  increase  the  altitude  of  the  peak  electron  concentration,  but  are  effective 
in  bringing  substantial  amounts  of  ionization  to  higher  altitudes . 

A quantitative  model  of  the  Spruce  ion  cloud  is  developed  that  is  in 

substantial  agreement  with  most  of  the  reliable  measurements  of  barium  ion 

clouds  produced  from  48  kg  chemical  payload  releases  at  around  190  kilometers 

altitude.  The  best  estimate  of  the  total  ion  inventory  developed  here  for  the 
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Spruce  ion  cloud  is  6 ± 1 x 10  . Detailed  models  A,  B and  C are  developed 

in  order  to  illustrate  the  variations  produced  by  different  assumptions.  The 

basic  model  A has  a scale  size  transverse  to  the  magnetic  field  of  2.66  km  and 
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a total  field-line  content  through  the  center  of  the  cloud  of  2.4  x 10  m . It 

13  -3 

reaches  its  maximum  electron  concentration  of  2.45  x 10  m at  approxi- 
mately 30  seconds  after  release.  The  subsequent  expansion  and  descent  of  the 
ion  cloud  is  shown  gr^hically.  Model  B has  the  same  initial  conditions  as 
model  A but  shows  that  downward  neutral  winds  and/or  a decreasing 
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neutral 'density  scale-height  can  maintain  the  peak  electron  concentration 
12  -3 

above  7 x 10  m for  longer  than  30  minutes  after  release.  Model  C 
presents  an  underestimate  of  the  total  ion  inventory  that  is  forced  to  agree 
precisely  with  two  early -time  radar  measurements  of  the  peak  electron  densi- 
ties which  in  model  A are  treated  as  lower  bounds. 
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1.  INTRODUCTION  AND  P’tmMARY 


In  order  to  model  the  characteristics  of  striations  that  have  developed 
from  barium  releases,  it  is  important  to  know  the  original  configuration  of 
ionization  in  the  ion  cloud  out  of  which  the  striations  are  formed.  The  electron 
densities  and  distribution  of  scale  sizes  that  exist  in  the  striations  at  late  time 
when  they  are  well -developed  depend  on  these  initial  conditions.  It  is  the  pur- 
pose of  this  report  to  describe  a procedure  for  developing  a model  of  the 
initial  ion  cloud  and  to  apply  that  procedure  to  the  Spruce  ion  cloud  that  was 
released  February  1,  1971  as  part  of  the  Secede  n Test  Series. 

As  described  in  more  detail  below,  three  aspects  of  developing  a 
description  of  the  ionization  in  the  ion  cloud  are  treated  in  the  next  three 
sections  of  this  report.  First,  the  expansion  of  the  barium  neutral  cloud 
which  consists  almost  entirely  of  vaporized  barium  atoms  is  treated.  Second, 
the  deposition  of  ionization  resulting  from  photoionization  of  the  barium  atoms 
in  the  expanding  neutral  cloud  is  calculated.  These  two  topics  are  treated  in 
Sections  2 and  3 and  result  in  a description  of  the  characteristics  of  the  ion 
cloud  at  early  time.  The  behavior  of  the  cloud  at  late  time,  that  is,  its  trans- 
port due  to  neutral  winds  and  electric  fields  and  its  change  in  shape  along  mag- 
netic field  lines  and  descent  into  the  neutral  atmosphere,  is  treated  in  Section 
4.  In  the  last  section  of  the  report  the  results  of  these  modeling  procedures 
are  applied  to  develop  quantitative  models  of  the  Spruce  ion  cloud. 

We  will  be  treating  several  aspects  of  the  continuity  equation  for  ions 
written  in  the  form 


V are  components  of  velocity 
z . 


■where  n is  the  ion  concentration,  v.  and 

U 

transverse  and  parallel  to  the  magnetic  field  and  S(r,  t)  is  the  source  function 
describing  the  creation  of  ionization.  In  Section  2 we  describe  a procedure  for 
calculating  the  source  function  S(r,  t)  . In  Section  3 we  neglect  the  second 
term  on  the  left-hand  side,  that  is,  the  transverse  transport  of  ionization 
across  magnetic  field  lines,  and  solve  the  remaining  parts  of  the  continuity 
equation  at  early  time  during  the  first  minute  following  release  when  S(r,  t) 
is  nonzero.  In  Section  4 we  describe  the  behavior  of  the  ion  cloud  beyond  one 
minute  after  release  when  the  right-hand  side  of  Eq.  (1. 1)  is  equal  to  zero 
under  the  assumption  that  the  ion  cloud  moves  as  a whole  across  magnetic 
field  lines. 


There  are  several  motivations  for  carrying  out  the  investigation 
described  in  this  report.  These  motivations  can  be  best  expressed  by  asking 
several  questions.  How  can  the  apparent  large  size  of  the  photographic  image 
of  neutral  clouds  be  reconciled  with  the  release  of  only  reasonable  amounts  of 
barium  vapor?  How  can  radar  measurements  of  the  transverse  dimension  of 
the  early  -time  ion  cloud  be  consistent  with  the  larger  size  of  the  expanding 
neutral  cloud?  How  quickly  can  ion  clouds  descend  to  low  altitudes  and  can  the 
peak  ion  concentration*  increase  at  late  time?  The  analyses  carried  out  in 
Sections  2,  3 and  4 provide  satisfactory  answers  to  these  questions. 

In  Section  2 a description  of  the  expanding  neutral  barium  atom  vapor 
cloud  is  provided.  The  expansion  proceeds  via  a snowplow  model  with  diffusion 
added.  The  inclusion  of  diffusion  is  a new  feature  and  provides  a description 
of  the  transition  of  the  initial  r£q>id  expansion  of  an  arbitrary  initial  density  con- 
figuration into  a Gaussian  profile  that  is  expanding  by  diffusing  through  the 
ambient  atmosphere  at  late  time.  The  effect  of  diffusion  is  to  change  an  arbi- 
trary initial  density  profile  into  a Gaussian  profile. 


We  have  carried  out  a detailed  analysis  of  the  intensity  profile  recorded 

on  a photographic  film  of  the  Spruce  neutral  cloud  at  20  seconds  after  release. 

♦Ellectron  concentration  and  ion  concentration  are  equal  to  each  other  and  will 
be  used  interchangeably  throughout  this  report.  Electron  density  or  ion  density 
is  sometimes  used  to  mean  the  same  quantity. 
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This  profile  is  compared  with  the  shapes  that  are  predicted  from  our 
theoretical  treatment  of  snowplow  plus  diffusion.  We  find  that  both  the  scale 
size  and  the  shape  predicted  by  the  theoretical  model  are  in  reasonable  agree- 
ment with  the  photographic  record  when  interpreted  properly.  Hence,  we 
believe  that  we  have  provided  a plausible  explanation  for  the  apparent  large 

photographic  image  of  neutral  clouds  during  the  first  tens  of  seconds  after 

1 2 

release.  Previous  theoretical  treatments  of  this  problem  ’ have  assumed 
that  unreasonably  large  amounts  of  neutral  barium  atom  vapor  were  produced 
in  order  to  match  the  size  of  the  early -time  images  of  neutral  clouds. 

Section  3 provides  a description  of  the  deposition  of  ionization  by  the 
photoionization  of  the  neutral  barium  atoms  in  the  expanding  neutral  cloud. 

As  a result  of  this  procedure  it  is  found  that  the  scale  size  of  the  ion  cloud 
transverse  to  the  magnetic  field  is  considerably  smaller  than  the  neutral  cloud 
radius.  The  central  high -density  portion  of  the  ion  cloud  is  more  peaked  than 
a Gaussian  distribution  but  the  outer  edges  fall  off  less  rapidly  than  do  a 
Gaussian  and  contain  a far  greater  inventory  of  ions. 

In  Section  4 several  aspects  of  the  motion  of  barium  ion  clouds  are 
treated.  General  equations  are  derived  for  the  motion  of  both  low-  and  high- 
conductivity  ion  clouds.  These  equations  incorporate  several  features  not 
ordinarily  included  that  have  a bearing  on  the  subject  of  the  behavior  and 
motion  of  the  ion  cloud.  They  include  a neutral  atmosphere  whose  density  is 
a function  of  altitude,  and  the  presence  of  a neutral  wind  that  is  not  necessar- 
ily spatially  uniform.  The  effect  of  gravity  and  a magnetic  field  with  an  arbi- 
trary dip  angle  is  also  included. 

Low-conductivity  ion  clouds  are  transported  across  magnetic  field 
lines  solely  by  ambient  electric  fields  which  produce  E x S motion.  The 
presence  of  a neutral  wind  tends  to  slide  the  ion  cloud  along  the  magnetic  field 
lines.  Neutral  winds  have  a greater  effect  on  the  transport  of  higher  conduc- 
tivity ion  clouds  across  magnetic  fields.  Indeed,  in  the  limit  of  very  high  con- 
ductivity, the  ion  cloud  moves  as  though  the  magnetic  and  electric  fields  were 
not  present.  The  magnetic  field  serves  solely  to  align  the  ion  cloud. 
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Both  gravity  and  the  component  of  il:e  neutral  wind  parallel  to  the 
magnetic  field  have  a large  effect  on  the  distribution  of  ionization  parallel  to 
the  magnetic  field.  The  changes  in  the  distribution  of  ionization  and  descent 
of  the  ion  cloud  are  computed  explicitly  for  various  examples  of  vertical  neutral 
winds.  It  is  found  that  downward  winds  decrease  the  rate  at  which  the  ion  cloud 
spreads  parallel  to  the  magnetic  field  and  thereby  lessens  the  decrease  in  the 
peak  ion  content. 

In  Section  5 we  apply  these  modeling  procedures  to  the  Spruce  ion 

cloud.  A quantitative  model  is  developed  that  is  in  substantial  agreement  with 

most  of  the  reliable  measurements  of  barium  ion  clouds  produced  from  48 

chemical  payload  weight  releases  at  around  190  km  altitude.  The  best  estimate 

24 

of  the  total  ion  inventory  developed  here  for  the  Spruce  ion  cloud  is  6 i 1 x 10  . 

Detailed  models  A,  B and  C are  developed  in  order  to  illustrate  the  variations 
produced  by  different  assumptions.  The  basic  model  A has  a scale  size  trans- 
verse to  the  magnetic  field  of  2.66  km  and  a total  field-line  content  through  the 

17  -2 

center  of  the  cloud  of  2.4  x 10  m .It  reaches  its  maximum  electron  con- 

13  -3 

centration  of  2.45  x 10  m at  approximately  30  seconds  after  release.  The 

subsequent  expansion  and  descent  of  the  ion  cloud  is  shown  graphically.  Model 

B has  the  same  initial  conditions  as  model  A but  shows  that  downward  neutral 

winds  and/or  a decreasing  neutral -density  scale -height  can  maintain  the  peak 
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electron  concentration  above  7 x 10  m for  longer  than  30  minutes  after 
release.  Model  C presents  an  underestimate  of  the  total  ion  inventory  that  is 
forced  to  agree  precisely  with  two  early -time  radar  measurements  of  the  peak 
electron  concentration  which  in  model  A are  treated  as  lower  bounds. 


2.  NEUTRAL  CLOUD  EXPANSION 


2.1  INTRODUCTION 

In  order  to  properly  model  the  characteristics  of  the  ion  cloud  (size, 
electron  concentration,*  etc)  produced  from  barium  releases,  it  is  necessary 
first  to  model  the  behavior  of  the  expanding  neutral  cloud.  The  distribution 
of  ionization  in  the  ion  cloud  can  then  be  calculated  by  using  the  neutral  cloud 
distribution  as  a source  for  the  production  of  barium  ions. 

In  an  earlier  work  ^ the  expansion  of  the  neutral  cloud  was  modeled 
in  three  separate  stages.  In  the  first  stage,  the  gas  release  is  slowed  down 
solely  by  losing  momentum  to  the  surrounding  air  via  a snowplow  model.  In 
the  second  stage  the  neutral  cloud  continues  to  expand  depending  on  the  range 
or  mean  free  path  of  the  barium  atoms  in  the  ambient  air.  After  a few  seconds 
during  which  the  initial  kinetic  energy  of  expansion  has  been  dissipated,  the 
neutral  atoms  continue  to  spread  by  diffusing  through  the  ambient  air. 

There  are  a number  of  fundamental  problems  associated  with  some  of 
the  assumptions  that  were  made  in  developing  that  model.  First,  it  was  as- 
sumed that  the  vaporization  efficiency  of  barium  atoms  was  30%  of  the  available 
excess  barium  contained  in  the  chemical  payload.  This  efficiency  is  greater 
than  has  ever  been  reliably  measured.  Furthermore,  it  was  assumed  that  the 
mass  of  the  expanding  neutral  cloud  was  20%  of  the  total  chemical  payload 
weight  amounting  to  a mass  j^proximately  40%  more  than  the  mass  of  the  bar- 
ium vapor.  There  is  no  other  species  in  significant  amounts  other  than  small 
amounts  of  the  impurity  strontium  to  account  for  this  extra  mass.  These 
unjustified  assumptions  were  motivated  by  the  ^parent  large  optical  size  of 
the  neutral  cloud  at  early  time . We  will  discuss  in  detail  how  more  reasonable 
assumptions  regarding  the  mass  of  vt^r  pjroduced  can  be  consistent  with 
optical  data. 

*l^ee  footnote  on  page  8. 
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The  simplified  snowplow  model  adopted  earlier  did  not  properly 
conserve  momentum  and  energy . We  will  use  a more  appropriate  snowplow 
model  here.  The  earlier  description  of  the  transition  from  the  initial  snow- 
plow regime  to  the  later  diffusion  regime  was  at  best  heuristic.  We  will  for- 
mally introduce  diffusion  into  the  snowplow  model  which  allows  assumed  initial 
density  distributions  to  evolve  into  a Gaussian  shape. 

The  results  of  the  former  modeling  were  to  produce  ion  clouds  that 
contained  a larger  ion  inventory  and  had  characteristic  transverse  dimensions 
larger  than  now  believed  to  be  reasonable.  The  principal  success  of  the  for- 
mer modeling  effort  which  speaks  in  its  favor  was  that  it  did  predict  the  cor- 
rect change  in  size  of  the  ion  cloud  and  the  peak  electron  densities  achieved  by 
a series  of  different  size  releases  conducted  at  differoit  altitudes.  The  var- 
iation in  chemical  payloads  ranged  over  two  and  one-half  orders  of  magnitude 
and  the  background  atmospheric  density  varied  by  a factor  of  twenty.  We 
shall  show  how  the  modifications  in  the  description  of  the  expansion  of  the 
neutral  cloud  mentioned  above  lead  to  neutral  and  ion  clouds  whose  physical 
properties  are  consistent  with  available  data. 

2. 2 SNOWPLOW  MODEL  WITH  DIFFUSION 

We  have  investigated  the  early -time  behavior  of  the  expanding  neutral 

3 

gas  cloud  by  using  an  extension  of  a snowplow  model  introduced  by  Stuart  in 
1965 . Much  of  our  development  of  this  model  was  done  in  connection  with  other 
work  relating  to  the  release  of  gas  clouds  at  satellite  orbital  velocities,  but 
the  specific  modification  of  considering  the  effect  of  high  release  velocity  is 
not  relevant  to  this  problem. 

In  this  model  all  lengths  scale  with  a characteristic  length  called  the 
equal-mass-radius,  a^  . This  length  is  the  radius  of  a spherical  volume  of 
ambient  atmosphere  with  mass  equal  to  the  mass  of  the  released  gas,  M^  . 

In  terms  of  the  ambient  atmospheric  density,  p , the  equal-mass -radius  is 

a 

defined  by  a = (3M  /4ir  p . The  central  assumptions  of  Stuart's  model 

w O A 

are  that  a spherical  cloud  of  gas,  consisting  of  released  gas  plus  swept  up 


ambient  atmosphere,  retains  a self -similar  or  shape -preserving  density 
profile,  that  the  velocity  of  the  gas  is  radial,  and  that  the  velocity  increases 
linearly  with  radius  from  the  center  of  the  cloud. 


This  model  appears  to  be  ^propriate  when  3 < a^A  < 30  where  X 

is  the  ambient  mean -free -path  between  collisions.  Since  a <*1-2  km  for 

o 

48  kg  barium  releases  at  around  190  km  altitude  and  X ~ 100-200  m,  this 
model  which  treats  the  released  vapor  and  ambient  atmosphere  as  being  well- 
mixed  but  still  subject  to  a fluid  treatment  is  appropriate.  At  higher  altitudes 
where  a^  < 3X  a kinetic  treatment  would  be  necessary.  At  lower  altitudes 
where  a^  > 30  X , it  is  difficult  to  consider  the  released  and  ambient  species 
as  being  well-mixed  and  other  gasdynamic  phenomena  neglected  in  the  snow- 
plow treatment,  such  as  shock  formation  and  wave  generation,  would  be  dom- 
inant. 

In  this  work  we  modify  the  Stuart  model  to  include  the  fact  that  at 
late  time  the  gas  cloud  can  be  modeled  by  a Gaussian  density  profile  diffusing 
into  the  ambient  atmosphere.  The  result  is  a model  which  initially  follows  a 
Stuart  snowplow  expansion  and  evolves  into  a Gaussian  diffusing  cloud. 
Obviously,  we  must  relax  one  of  the  central  assumptions  of  the  Stuart  model, 
specifically  that  the  cloud  retains  a shape -preserving  density  profile.  We 
will  treat  the  effect  of  diffusion  as  if  it  were  superimposed  on  the  e3q)ansion 
process  modeled  by  the  self -similar  snowplow  expansion.  In  a self -similar 
expansion  the  size  of  the  cloud  is  given  by  R(t)  , a time -dependent  scale  length 
of  the  cloud.  If  the  self -similar  sh^  were  a spherical  constant -density  pro- 
file, for  example,  we  would  take  R(t)  to  be  equal  to  the  radius  of  the  sphere. 

We  use  this  same  R(t)  to  scale  the  continuity  equation  for  mass,  as 
was  used  in  the  development  of  Stuart's  model,  with  the  addition  of  a term  due 
to  diffusion; 


DP 

Dt 


+ p7 


= 


Mp 

M ^ 


(2.1) 


where  the  hydrodynamic  derivative  is  given  by 


I 
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D 

DT 


+ V • V , 


(2.2) 


' (•^)r 

V is  the  fluid  velocity  of  the  gas  given  by  ^ = rR/R  , where  R(t)  is  the 

characteristic  scale  length  of  the  expanding  cloud  if  there  were  no  diffusion, 

p is  the  mass  density  and  M(t)  is  the  total  mass  of  the  gas  cloud  including  the 

swept-up  ambient  mass.  The  dot  indicates  differentiation  with  respect  to  time. 

D is  the  diffusion  coefficient  which  was  not  included  in  earlier  treatments  of 
n 

the  Stuart  model.  The  last  term  on  the  right-hand-side  of  Eq.  (2.1)  includes 
the  effect  of  sweepii^  up  ambient  mass  according  to  Stuart's  original  model. 

By  making  a suitable  change  of  space  and  time  variables  we  can  define 
a function  f(^(r,  t),  rQ(t))  which  is  analogous  to  the  shape  function  of  the  self- 
similar  model.  If  we  define  4 and  by 


I = r/Rd) 


4D^dt 

R^(t) 


(2.3) 


implying  that  diffusion  begins  at  t = t^^  , then  the  density  can  be  written 

p(?,t)  = [M(t)/R^t)]f[?(?,t).  Tj5(t)]  . 

With  these  substitutions  Eq.  (2. 1)  can  be  simply  written 


(2.4) 


where  we  have  used  the  following; 


-»  R ■* 

V = ^r  ; 

V • V = 3 R/R  ; 
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Equation  (2.4)  is  the  diffusion  equation  for  f(^,  r)  , which  has  the 


solution 


f I.Tj,  = / f(|^0)— A ^ . 


(2.5) 


In  terms  of  this  function  we  may  write  the  density  at  any  time  t > t^  as 

,(?,  t)  = ^ fih'p  (r',  t„)  — 

where  Rjj  = R(tjj)  and  = M(tjj)  . For  « 1 we  have  the  limit  of 

little  diffusion  and  the  density  reduces  to  that  given  by  the  self -similar  model, 


or 


p(r,t)  = -^f(r/R,  0) 
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At  long  times  after  diffusion  has  been  turned  on  the  similarity  function  f 
approaches  a Gaussian  distribution  with  an  expanding  radius  and  constant 
volume: 

as  Tjj  00  . 

In  order  to  illustrate  the  effect  that  diffusion  in  i has  on  an  assumed 
density  distribution,  we  have  applied  Eq.  (2.5)  to  two  specific  examples.  The 
first  is  a uniform  density  out  to  a finite  expanding  radius;  we  refer  to  this  shape 
as  a flat -top  distribution.  The  second  is  an  expandii^  spherical  shell  with  all 
the  mass  concentrated  at  the  surface  of  the  shell.  For  the  flat -top  distribution 


f(€,  0)  = 3/4ff  ; 4 < 1 
f(«,  0)  = 0 ; 4 > 1 , 

and  for  the  spherical  shell 

f(C,  0)  = 6(€-l)/4ir  . 

The  normalization  has  been  chosen  so  that  each  contains  a unit  volume . 

We  have  found  that  the  approach  to  a Gaussian  profile  takes  place 

more  r^idly  if  the  Gaussian  radius  has  the  form  |o(  + instead  of 

T ^ L e 
D > > 

In  order  to  illustrate  the  change  of  the  flat -top  and  spherical-shell  distribution 
toward  the  Gaussian  form  given  in  Eq.  (2. 7)  at  various  values  of  , we  have 
renormalized  both  the  radius  and  density  by  defining 


z = {/ro.693(oi+ Tj,)]* 
g(z,  Tn)  = [z(»*rp)]^  f(«,  Tp) 


The  limiting  value  of  g(z,  t^)  for  large  Tjj  is  exp(-  0.693  . The  scale 

factor  0.693  is  chosen  so  that  g(l,  <»)  is  0.5. 

Figures  2. 1 and  2.2  show  g(z,  T^)  for  various  values  of  Tjj  for  the 
flat -top  and  spherical  shell  initial  distributions  respectively.  The  value  of  a 
in  the  scalings  given  by  Eq.  (2.8)  are  different  for  the  two  cases.  It  is  appro- 
priate that  a be  unequal  to  zero  because  in  both  cases  the  initial  shapes  have 
finite  spatial  extent.  We  have  chosen  et  so  that  g(l,  1)  is  approximately 
equal  to  one  half;  for  the  flat -top  at  = 0.5  and  for  the  spherical  shell  at  = 0.75. 
The  larger  value  of  at  for  the  spherical  shell  results  from  the  fact  that  all  of 
the  mass  is  initially  at  C = 1 for  the  shell  while  the  mass  is  distributed  over 
€ ^ 1 for  the  flat -top. 

The  vertical  dashed  lines  in  Figures  2. 1 and  2.2  indicate  the  position 
of  the  initial  radius  in  the  newly-scaled  length  coordinate  z . For  greater 
than  one,  the  position  of  the  initial  radius  moves  to  the  left  and  approaches 
z = 1.2  as  Tjj  indicating  the  increase  in  size  of  the  cloud  due  to 
diffusion.  For  the  flat -top  in  Figure  2.1,  the  horizontal  dashed  lines  indicate 
the  height  of  the  original  flat -top  in  the  newly-scaled  density  with  amplitude 
[ir(ai+  Tjj)]^f  («,  0). 

Examination  of  Figures  2. 1 and  2.2  indicates  that  the  density  distri- 
butions are  approaching  Gaussian-like  shapes  for  > 0.  2 and  0.  75  for  the 
two  cases,  respectiveiy,  and  have  essentially  the  asymptotic  shape  given  in 
Eq.  (2.7)  for  >0.4  and  1.0  respectively.  We  will  discuss  specific  eval- 
uations of  in  the  next  subsection. 

The  ^proximation  of  f(€,  by  its  limiting  form  [»(«+ 
exp  [-4  /(«  + fj^)]  is  equivalent  to  approximating  the  density  distribution 
of  the  cloud  by  a Gaussian  density  profile 
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The  evolution  of  the  normalized  density  distribution  from  an  initial  flat-top  to  ; 
Gaussian  as  a result  of  diffusion  for  several  values  of  the  diffusion  parameter 


01 


(2.9) 


p(r,t)  - M(t)  exp  [-  r^/a^(t)] 

vrith  a time -dependent  Gaussian  radius,  a(t),  given  by 

a^t)  = rV)  [«  + TjjCt)]  (2.10) 

where  is  given  by  Eq.  (2. 3).  Although  the  actual  shsqpe  at  early  time  is  not 
Gaussian,  in  the  case  of  the  flat -top  distribution  we  find  that  the  distance  at  which 
the  density  is  l/e  of  its  value  at  r = 0 is  within  a few  percent  of  R(t)  for  Tj^(t)  <0.  5 
and  of  course  is  indistinguishable  from  a(t)  given  by  Eq.  (2. 10)  for  > 0.  5. 

In  the  preceding,  we  have  discussed  a procedure  for  introducing 
diffusion  within  the  framework  of  a snowplow  model.  The  effect  of  diffusion 
is  to  transform  an  initial  self-similar  profile  for  the  density  into  a Gaussian 
profile.  The  treatment  so  far  does  not  depend  on  the  time-dependence  of  R(t) 
which  we  take  to  be  determined  by  the  snowplow  model  without  diffusion. 

Using  local  conservation  of  mass  and  momentum  and  global  conservation 
of  energy,  the  Stuart  snowplow  determines  a function  0(s)  such  that 

R(t)  = 3^0(8)  (2.11) 

and 

s = Rot/a^  (2.12) 

where  R^  is  the  initial  expansion  velocity  of  R(t)  . 

Our  normalizations  for  and  s are  different  from  those  of  Stuart; 
with  our  normalizations  0 satisfies  the  differential  equation; 

where  y is  the  ratio  of  specific  heats.  The  solution  of  this  equation  has  the 
following  properties  for  large  and  small  s : 
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0(s)  -»  s for  s <<  1 

/ i\°-2  n 4 

0(s)  -♦  [6.25 8 for  s » 1 

At  early  time  the  cloud  expands  with  the  initial  expansion  velocity  R while 

0. 4 ” 

at  later  time  the  cloud  radius  is  predicted  to  expand  as  t ‘ .At  very  late 
time  when  the  snowplow  description  is  no  longer  appropriate,  the  cloud 
expands  as  t®'®  due  to  diffusion.  It  is  apparent  that  once  the  appropriate 
value  of  y of  the  released  gas  is  specified,  the  snowplow  model  for  the 
expansion  of  a gas  cloud  determines  the  time -dependence  of  the  radius  R{t) 
of  the  cloud  in  terms  of  the  two  key  parameters  a^  and  R^  . 

Once  the  functional  dependence  of  <p{s)  is  known,  then  we  may  define 
a normalized  diffusion  function  f (s)  by 

f(s)  = f . (2.13) 


The  diffusion  parameter  defined  in  Eq.  (2.3)  can  be  expressed  in  terms 
of  f(s)  as 

4D 

’■d<«  * 

where  s = R t/a  and  the  diffusion  is  assumed  to  commence  at  t = L . 

o o w 

For  a gas  with  y = 5/3  which  is  ^propriate  for  barium  atom  vapor,  we  show 
the  resulting  normalized  diffusion  function  f(s)  in  Figure  2.3.  Tht  solid 
curve  is  the  correct  value  of  f(s)  only  as  long  as  0(s)  is  given  by  the  snow- 
plow model.  If  the  snowplow  model  were  to  be  turned  off  at  a specific  value  of 

O 

s , then  0 in  the  denominator  would  be  a constant  and  <f>  would  increase 
linearly  with  s . Examples  are  shown  by  the  dashed  curves  in  Figure  2.3 
if  the  snowplow  were  stopped  at  s = 18  and  46  . We  will  discuss  specific 
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Figure  2.3 
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s - K„  t/a„ 
o o 

Normalized  diffusion  function  f(s)  for  a gas 
with  y = 5 /3  in  terms  of  which  the  diffusion 
parameter  is  defined  according  to  Eq.  (2.14). 
The  dashed  curves  show  the  values  of  f(s)  if  the 
snowplow  were  stopped  at  s = 18  and  46  respec- 
tively . 
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numerical  values  for  the  expansion  and  diffusion  of  the  Spruce  neutral  cloud 
in  the  following  subsection. 

2.3  EXPANSION  AND  SHAPE  OF  THE  NEUTRAL  CLOUD 

In  Section  5 of  this  report  we  will  present  parameters  that  define  two 
early -time  models  of  the  Spruce  ion  cloud.  Model  A and  model  C represent 
ion  clouds  whose  ion  inventories  are  an  upper  estimate  and  an  underestimate, 
respectively . The  model  parameters  listed  there  are  believed  to  be  valid 
because  they  are  consistent  with  the  model  of  the  neutral  gas  cloud  expansion 
described  in  the  previous  subsection  and  consistent  with  a detailed  analysis 
of  a film  record  of  the  Spruce  neutral  cloud  to  be  described  here. 

In  order  to  compare  the  results  of  the  theoretical  description  with 
the  observation  of  the  Spruce  neutral  cloud,  we  will  use  the  values  appropriate 
for  model  A.  A discussion  of  the  sensitivity  of  the  adopted  model  values  to 
various  assumptions  will  be  deferred  to  Section  5 . 

The  Spruce  cloud  was  released  on  1 February  1971  at  2352  UT  at  an 
altitude  of  about  190  km  from  Eglin  AFB  in  Florida.  The  CIRA  1965  atmos- 
pheric model  6 at  1800  hours,  appropriate  for  the  solar  flux  conditions  that 
existed  durii^  the  second  half  of  January  1971,  was  used  to  obtain  the  follow- 
ing parameters  at  190  km  altitude: 

-3 

ambient  mass  density  p = 0.45  kg  Am  ; 

A 

neutral  air  temperature  T^^  = 1086*  K ; 

barium  ion -neutral  collision  time  t = 0.8  s ; 

2 

barium -neutral  diffusion  coefficient  D^  = k T^  t/M  = 0.052  km  /s 
The  barium-neutral  collision  time  is  discussed  in  Section  4.4. 

The  chemical  formulation  for  the  barium  thermite  is  based  on  the 
reaction  of  2 . 5 moles  of  barium  per  mole  of  cupric  oxide  according  to  the 
formula 

2.5  Ba  + CuO  -♦  BaO  + Cu  + 1.5  Ba  . 
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In  addition,  1.8%  of  the  thermite  weight  was  barium  azide.  ^ As  a result, 

roughly  half  (47.8%)  of  the  chemical  payload  weight  was  available  as  excess 

26 

barium,  amounting  to  1.01  x 10  barium  atoms  per  48-kg  barium  release. 

The  results  of  model  A for  the  ion  cloud  discussed  in  Section  5 indicate  that 
approximately  15.5%  of  the  available  barium,  or  = 3.56  kg,  exits  the 
caimister  as  barium  atoms . We  note  that  this  corresponds  to  an  efficiency  by 
weight  of  7.4%  which  is  right  in  the  middle  of  the  usual  estimate  of  5 to  10%. 
Measurements  of  the  initial  expansion  velocity  report  typical  values  between 
1.0  and  1.4  km/s.  For  model  A we  have  adopted  the  value  1. 15  km/s  for  the 
initial  expansion  velocity  . 

We  have  now  specified  the  parameters  necessary  to  apply  the 

theoretical  model  developed  in  the  previous  subsection  to  observations  of  the 

neutral  cloud.  By  using  the  above  mass  of  released  vapor,  we  determine  the 

equal-mass-radius  for  model  A to  have  the  value  a^  = 1.  26  km.  The  dimen- 

o 

sionless  coefficient  of  the  diffusion  parameter  in  Eq.  (2. 14)  has  the  numer- 
ical value  4D^/a„R„  = 0.144. 
n o o 

The  solid  curve  in  Figure  2.4  shows  the  radius  of  the  Spruce  neutral 
cloud  as  a function  of  time  after  release  according  to  the  snowplow  model. 

Since  the  radius  and  time  are  given  on  a logarithmic  scale,  the  solid  curve  also 
represents  the  function  0(s)  shifted  horizontally  by  the  scale  factor 
and  vertically  by  the  factor  a^  . We  have  found  that  the  function 

9J(s)  = (l.44  s°-®  - 0.63)*  (2.15) 

is  indistinguishable  from  <p{s)  for  s > 1 and  have  used  0(s)  in  carrying  out 
calculations  involving  <^(s)  . The  radius  of  the  neutral  cloud  as  a function  of 
time  is  then  given  by 


->* 


0.63 


(2. 16) 
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Neutral  gas  cloud  r^us  as  a function  of  time  after  release  according  to  model  A with 
ju  = 1.26  km  and  = 1.15  km/s.  The  dashed  curves  show  the  effect  of  suddenly  turning 
da  the  snowplow  at  20  s and  50  s. 


The  neutral  cloud  reaches  the  equal -mass -radius  of  1.26  km  at 
s = 1.16  or  at  around  1.3  seconds  after  release.  At  this  time  the  expanding 
barium -vapor  cloud  has  picked  up  a mass  of  ambient  air  equal  to  the  mass  of 
the  expanding  vsq>or.  It  'would  seem  inappropriate  to  introduce  diffusion  as  an 
important  mechanism  for  times  earlier  than  this  time.  For  illustrative  pur- 
poses we  have  started  diffusion  a little  later  at  around  2 seconds  after  release 
when  s has  the  value  of  1.8. 

According  to  the  discussion  in  the  previous  subsection,  diffusion 
initially  does  not  increase  the  size  of  the  neutral  cloud.  Until  = 0.5  , 
it  serves  primarily  to  evolve  an  initial  flat -top  distribution  into  a Gaussian 
shape.  During  this  time  the  radius  of  the  cloud  given  by  the  solid  curve  or 
Eq.  (2. 16)  gives  the  dimension  of  the  point  on  the  neutral  cloud  that  has  a value 
1/e  of  its  value  at  r = 0. 

The  dashed  curves  in  Figure  2.4  show  the  effect  of  suddenly  turning 
off  the  snowplow  model  at  t = 20  s and  t = 50  s after  release.  According 
to  the  snowplow  model  the  neutral  cloud  1/e  radius  has  a value  of  4. 72  km  at 
20  s and  the  dimensionless  time  parameter  s has  the  value  18  . The  diffu- 
sion parameter  can  be  evaluated  by  using  its  definition  in  Eq.  (2. 14),  the 
numerical  value  0.144  for  the  coefficient  4 Dj^/a^  , and  Figure  2.3  for  f(s). 
We  find  that  is  approximately  0. 4 at  20  s after  release  and  reaches  the 
value  of  0.5  at  around  35  s after  release.  We  should  remark  that  the  solid 
curve  in  Figure  2.3  is  valid  only  during  the  time  that  the  snowplow  is  describ- 
ing the  time  dependence  of  R(t)  . 

If  the  snowplow  were  turned  off  at  20  s after  release,  diffusion  would 
cause  the  neutral  cloud  to  continue  to  evolve  toward  a Gaussian  but  would  not 
increase  its  size  for  the  next  10  to  15  seconds.  At  this  time  when  reached 
the  value  0.5,  the  neutral  cloud  'would  begin  to  increase  in  size.  On  the  other 
hand,  if  we  delayed  the  turnoff  of  the  snowplow  until  50  s after  release,  diffu- 
sion would  increase  the  size  of  the  neutral  cloud  between  35  and  50  s as  shown. 
After  the  sudden  turnoff  at  50  s , the  neutral  cloud  would  continue  to  expand  by 


diffusion.  We  note  that  both  dashed  curves  are  approaching  the  dashed  curve 
labeled  asymptote  which  at  late  time  gives  the  Gaussian  radius  of  the  neutral 
cloud  as  4D^t  . Of  course,  the  sudden  turnoff  of  the  snowplow  is  unphysical; 
in  actuality  it  gradually  ceases  to  contribute  to  the  increase  in  size  of  the 
neutral  cloud.  Measurements  of  the  images  of  neutral  clouds  on  photographic 
plates  tend  to  show  that  during  the  time  period  from  20  to  100  seconds  after 
release  they  do  not  increase  in  size  very  rapidly.  We  believe  that  this  is  a 
reflection  of  the  fact  that  the  snowplow  model  is  no  longer  an  appropriate 
description  for  the  increase  in  size  of  the  neutral  cloud  at  this  late  time. 

In  Section  5,  we  will  discuss  the  increase  in  size  of  the  ion  cloud. 

For  reference  purposes  the  dotted  line  in  Figure  2.4  shows  the  Gaussian  half- 
length  of  the  ion  cloud  parallel  to  the  magnetic  field.  This  curve  is  derived 
by  assuming  that  the  ion  cloud  increases  in  length  according  to  the  formula 

t^t)  = (4.72)^  + 4Djj(t-20)  km^  (2.17) 

where  is  the  parallel  diffusion  coefficient, 

Do  - k(To  * T,)r/M  . ((T^  + T,)/T„)D„  . 

For  ion  temperatures  equal  to  the  neutral  temperature  and  electron  tempera- 
tures 30%  hotter,  D has  the  value  0. 12  km^/s  . 
o 

2.4  COMPARISON  WITH  OBSERVATION  OF  THE  SPRUCE  NEUTRAL 

CLOUD 

0 

Boquist  et  al.  have  reported  a measurement  of  the  radius  of  the 
Spruce  neutral  cloud  at  20  seconds  after  release  of  7.2  km.  This  measured 
size  is  shown  as  a dot  in  Figure  2.4.  This  measurement  is  seen  to  be  consid- 
erably larger  than  the  1/e  value  predicted  by  the  snowplow  plus  diffusion  model. 
This  result  is  typical.  Generally,  the  measured  diameter  of  the  well-exposed 
photographic  image  of  neutral  clouds  are  considerably  larger  than  would  be 
predicted  by  the  snowplow  model. 


27 


As  we  shall  see  in  the  next  section,  the  size  of  the  neutral  cloud  at  early 
time  has  a considerable  impact  on  the  deposition  of  ions  in  the  ion  cloud.  We 
therefore  undertook  a detailed  examination  of  the  photographic  image  of  the 
Spruce  neutral  cloud  at  20  seconds  after  release.  Technology  International 
Corporation  kindly  supplied  us  with  a densitometer  tracing  of  a horizontal  scan 
of  the  Spruce  neutral  cloud  at  20  seconds  after  release  from  film  record  #71715. 
This  film  is  a color  positive  film  so  that  the  image  of  the  neutral  cloud  appears 
light-colored  compared  to  the  dark  blue  sky  background. 

Figure  2.  5 is  a representation  of  the  measured  transmittance  of  the 
image  of  the  neutral  cloud.  The  image  is  slightly  asymmetric  and  the  dashed 
vertical  line  is  drawn  through  the  peak  transmittance.  Also  indicated  by  the 
short  vertical  bars  and  the  horizontal  arrow  is  a quantitative  representation  of 
the  full  width  of  the  Spruce  neutral  cloud  that  had  been  previously  measured  by 
eye.  We  see  that  this  measured  width  corresponds  to  the  outer  fringes  of  the 
neutral  cloud  where  it  begins  to  have  a significant  contrast  with  respect  to  the 
background  sky.  One  millimeter  on  the  film  plane  corresponds  to  1. 16  km. 

The  characteristics  of  exposed  film  are  usually  given  by  D-log  E curves 
where  the  optical  density  D = log  (1/T),  T is  the  transmittance,  and  E is  the 
energy  density  (erg/cm  ) incident  on  the  film.  Such  characteristics  are  not 
available  for  color  film  in  general.  We  have  attempted  to  extract  useful  quanti- 
tative information  regarding  the  actual  distribution  of  light  intensity  incident  on 
the  film  by  trying  various  models  of  film  behavior  that  produce  reasonable  D- 
log  E curves.  These  various  analyses  produced  similar  qualitative  results  and 
the  quantitative  differences  were  not  great.  As  an  example,  if  the  number  of 
absorbers  remaining  in  the  exposed  film  decreases  exponentially  with  the  number 
of  photons  incident,  then  the  density  of  the  film  after  development  is  given  by 
D = Dp  ejq)  ( - k I)  where  k and  are  arbitrary  constants.  The  light  intensity 
incident  on  the  film  can  be  written  I = + Ip(r)  where  is  the  intensity  of 

the  sky  background  and  I (r)  is  the  intensity  of  light  from  the  cloud.  The  dis- 
tribution  of  intensity  of  the  neutral  cloud  can  then  be  expressed  in  terms  of  the 
transmittance  of  the  image  of  the  neutral  cloud  as 
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Figure  2.5.  Measured  transmittance  of  the  Spruce  neutral  cloud  image  at 
20  seconds  after  release  as  recorded  on  film  record  #71715. 
(Original  densitometer  tracing  courtesy  of  W.  Boquist, 
Technology  International  Corporation.) 
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(2. 18) 


yr)  = Clog 


log(l/T) 

“>e(‘^b) 


where  is  the  transmittance  of  film  background. 


The  original  of  Figure  2. 5 was  used  to  determine  the  intensity  of  the 
image  of  the  neutral  cloud  as  a function  of  radius  according  to  Eq.  (2. 18).  Be- 
cause of  the  asymmetry  in  the  image  of  the  neutral  cloud,  three  different  distri- 
butions were  determined.  One  distribution  corresponded  to  the  right -half  of  the 
cloud,  the  second  corresponded  to  the  left -half  of  the  cloud  and  the  third  corres- 
ponded to  the  half -width.  For  the  right-half  and  left -half  of  the  cloud,  the  appro- 
priate measured  background  transmittances,  0. 18  and  0. 19  respectively,  were 
used.  The  third  determination  provided  a different  distribution  than  the  average 
of  the  first  two  because  the  right-half  value  0. 18  was  used  for  the  background 
transmittance.  The  D-log  E film  characteristic  that  results  in  Eq.  (2. 18)  is 

consistent  with  a ratio  of  I (0)/l|,  = 2.  5 which  is  in  agreement  with  typical 

® ° -2  -1  -1  ** 
measured  values  of  I (0)  ~ 0. 3 erg  cm  s sr  and 

s"*  sr“  at  a solar  depression  angle  near  8 . 


Ijjj  ~ 0. 125  erg  cm 


Figure  2, 6 shows  a plot  of  the  logarithm  of  the  intensity  in  arbitrary  units 
as  determined  above  versus  the  square  of  the  radius  for  the  three  different  sets 
of  measurements.  If  the  transmittance  of  the  film  were  related  to  the  intensity 
in  the  cloud  by  Eq.  (2. 18),  and  if  the  intensity  distribution  in  the  cloud  were 
Gaussian,  the  data  points  would  lie  on  a straight  line.  We  see  that  for  more 
than  an  order  of  magnitude  change  in  intensity  beginning  at  approximately  one- 
half  the  value  of  the  central  intensity,  the  distribution  does  appear  to  be  Gaussian. 


We  have  attempted  to  compare  the  scale  length  and  the  shape  determined 
by  these  data  points  with  the  same  quantities  determined  by  the  theoretical  model 
developed  earlier  in  this  section.  With  the  assumption  that  the  l/e  point  is  given 
by  Eq.  (2. 16)  and  has  reached  a value  of  4. 7 km  at  20  s after  release,  and  that 
diffusion  was  turned  on  at  jyiproximately  2 s after  release,  we  have  determined 
the  resulting  shapes  of  the  model  density  distribution  for  values  of  around 
0. 4 ^ropriate  for  20  s after  release.  Three  such  sh^es  determined  by  the 
theoretical  model  are  shown  by  the  three  solid  curves.  The  central  intensity  of 
each  curve  was  adjusted  arbitrarily  in  order  to  match  the  data  points  at  large 
radius. 
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INTENSITY,  I (r),  (ARBITRARY  UNITS) 
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Figure  2.6.  Comparison  of  the  measured  intensity  distribution  of  the  Spruce 
neutral  cloud  at  20  seconds  after  release  with  the  electron 
density  profile  determined  from  a snowplow  model  with 
diffusion  added. 
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The  comparison  of  the  theoretical  model  shapes  with  the  shs^es  determined 
from  the  photogr^h  indicates  excellent  agreement  for  values  of  intensity  roughly 
a factor  of  2 below  the  measured  peak  intensity  of  the  cloud.  Two  important  fea- 
tures of  the  theoretical  curves  appear  to  be  consistent  with  the  data.  The  outer 
edges  of  the  density  profile  are  fit  well  by  a Gaussian  shape,  exp  (-  r^A>^  , and 
the  values  of  the  Gaussian  scale  length,  b,  (indicated  in  Figure  2. 6)  are  in  agree- 
ment. This  agreement  indicates  that  a value  of  the  diffusion  parameter  of 
around  0. 4 is  z^propriate.  We  further  note  that  the  scale  length  4.  7 km  deter- 
mined by  the  snowplow  model  as  the  radius  at  which  the  density  has  fallen  to  1/e 
of  its  central  value  is  consistent  with  the  observed  intensity  distributions. 

We  feel  that  this  analysis  provides  a plausible  explanation  for  the  apparent 
discrepancy  in  Figure  2. 4 in  comparing  the  measured  cloud  radius  with  the  results 
of  the  theoretical  modeling  according  to  snowplow  with  diffusion  added.  It  is  sug- 
gested that  the  dimensions  of  well -exposed  images  of  neutral  clouds  during  the 
first  few  minutes  after  release  (and  after  the  first  few  seconds)  greatly  exceed 
the  Gaussian  radii  of  the  clouds. 

The  measured  intensity  at  the  center  of  the  ion  cloud  is  somewhat  depressed 
below  the  theoretical  model  values.  There  are  three  possible  factors  contributing 
to  this  observation.  One  possibility  is  that  the  central  portion  of  the  neutral  cloud 
was  less  peaked  at  20  s after  release  than  predicted  by  the  theoretical  model. 

This  flatter  shape  would  be  consistent  with  the  assumption  that  the  emitted  dis- 
tribution of  barium  vapor  at  these  low  altitudes  is  more  shell-like  than  the  flat- 
top  distribution  we  have  assumed.  It  is  more  likely  that  the  derived  intensity 
distribution  is  not  proportional  to  the  actual  density  distribution  within  the  cloud. 
First,  c^tical  depth  effects  within  the  cloud  would  lead  to  a flatter  intensity  dis- 
tribution than  the  actual  distribution  of  vapor.  Second,  Eq.  (2. 18)  may  not  be 
reasonable  in  the  highly  exposed  regions  where  the  film  may  be  approaching  sat- 
uration. Indeed,  a simple  modeling  of  the  saturation  effect  results  in  a more 
peaked  distribution  of  the  data  points  in  the  central  region  without  changing  the 
data  points  in  the  outer  region  significantly.  The  corresponding  best  fit  value  of 
would  be  reduced  to  0. 3-0. 35  which  is  still  within  a reasonable  range. 


3.  DEPOSITION  OF  THE  ION  CLOUD 


In  this  section  we  discuss  the  creation  of  the  ion  cloud  from  the 
expandii^  neutral  barium  atom  cloud  that  was  discussed  in  the  previous 
section.  We  will  describe  the  procedure  for  calculating  the  ion  deposition 
and  defer  the  development  of  a specific  model  for  Spruce  until  Section  5 . 

In  this  section  we  are  concerned  only  with  the  early  development  of 
the  ion  cloud,  that  is,  with  the  deposition  of  the  ionization  which  takes  place 
during  the  first  minute  following  release.  In  Section  4 we  will  describe  the 
motion  of  the  ionization  and  its  change  in  shape  along  magnetic  field  lines  at 
times  greater  than  one  minute  after  release.  In  Section  5 we  will  put  together 
a unified  model  for  the  Spruce  ion  cloud  by  combining  the  procedures  developed 
in  this  section  for  the  deposition  of  the  ions  at  early  time  and  the  later  time 
behavior  of  the  ion  cloud  discussed  in  the  next  section. 

We  first  give  a qualitative  description  of  the  deposition  process.  By 
about  one  second  after  release  the  neutral  barium  atom  vapor  has  expanded  at 
essentially  constant  velocity  and  approaches  a spherical  volume  with  approx- 
imately an  equal-mass -radius.  At  this  time  the  concentration  of  the  barium 
atoms  within  the  neutral  cloud  is  comparable  to  the  ambient  local  concentration 
of  around  10  particles  m~  . Solar  radiation  excites  ground-state  barium 
atoms  into  metastable  states.  The  population  of  neutral  barium  atoms  is 
depleted  by  two  competing  processes:  photoionization  from  the  excited  meta- 

7 

stable  states  and  oxidation  by  ambient  oxygen  molecules.  Haerendel  of  the 
Max-Planck  Institute  has  consistently  reported  photoionization  times  of  order 

Q 

19  seconds  while  Best  of  AFCRL  reports  a time  closer  to  30  seconds.  In  our 
model  we  will  choose  24  seconds  as  being  within  20%  of  the  correct  value.  We 
shall  also  take  the  oxidation  rate  constant  of  barium  by  atmospheric  0«  to  be 
7.5  X 10  cm  /s  . The  fraction,  f , of  available  barium  atoms  that  are 
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ionized  is  equal  to  the  ratio  of  the  photoionization  rate  to  the  total  loss  rate 
of  barium  atoms  due  to  both  photoionization  and  oxidation, 

f = g . (3.1) 

1 + 1.8  X 10"®  N,- 
^2 


The  time  constant,  t.  , for  the  loss  of  Ba  is 

1 

T.  = 24f  = 5 s . (3.2) 

^ 1 + 1.8  X 10‘‘'n_ 

^2 

For  the  Spruce  cloud  released  at  190  km,  the  ionization  time  constant,  , 
is  equal  to  10  seconds. 

In  the  previous  section  we  discussed  the  expansion  of  the  neutral 
cloud.  We  have  shown  that  by  including  diffusion  the  neutral  cloud  quickly 
evolves  towards  a Gaussian  shape  with  an  equivalent  Gaussian  radius,  R(t)  , 
given  by  Eq.  (2. 16).  We  designate  by  N,j,  the  total  number  of  barium  atoms 
produced  in  the  release.  Since  only  the  fraction  f of  these  atoms  become  ions, 
and  the  time  scale  on  which  the  neutrals  are  depleted  is  , we  find  that  we 
can  write  the  source  function  for  the  creation  of  barium  ions  as 

S(r,t)  . ta(r,t)/r,  = ^ . (3.3) 

where  n^(r,  t)  is  the  concentration  of  barium  atoms. 

We  designate  the  field  line  content  yndz  by  N(r,t)  where  r now 
represents  a radius  perpendicular  to  S,  and  neglect  the  transport  transverse 
to  the  magnetic  field  during  the  ion  deposition.  By  integrating  the  continuity 
equation  (1.1)  along  magnetic  field  lines,  we  can  derive  the  following  equation 
which  governs  the  time  and  spatial  dependence  transverse  to  the  magnetic  field 
of  the  deposition  of  ions  on  a field  line; 
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It 


(3.4) 


expizrWw)  e-(t-tj)/T.  ^ 

R^(t) 


If  we  start  the  ionization  process  at  time,  t. , after  release  we  can  integrate  Eq. 
(3. 4)  in  order  to  determine  the  radial  profile  of  the  field-line-content  of  the  ion 
cloud; 


It  is  appropriate  that  we  delay  the  start  of  photoionization  until 
times  greater  than  1-1.5  seconds  when  the  neutral  cloud  has  expanded  to  an 
equal-mass -radius.  Prior  to  this  time  the  concentration  of  barium  atoms  is 
so  high  that  optical  depth  effects  prohibit  the  metastable  exciting  radi;i  ion  to 
permeate  the  neutral  cloud.  Indeed,  for  very  large  barium  releases  uniform 
excitation  of  the  neutral  barium  atoms  might  be  delayed  for  several  seconds. 
The  barium  atom  and  neutral  air  densities  only  approach  the  nominal  ambient 
values  at  times  longer  than  it  takes  R(t)  to  exceed  the  equal-mass-radius. 

On  the  other  hand,  radar  sensors  have  shown  that  the  peak  ion  concentration 

9 

rises  quite  rapidly  within  the  first  few  seconds  after  release  so  t^  cannot  be 
taken  to  be  too  large. 

To  exemplify  the  distribution  of  ionization  on  field  lines  that  results 
from  this  model  we  have  carried  out  the  necessary  integrations  to  determine 
the  field -line  ccxitent  as  a function  of  radius  transverse  to  the  magnetic  field. 
Explicitly,  if  we  define  a as  the  size  of  the  Gaussian  radius  of  the  neutral 
cloud  at  the  time  t = t.  + after  release,  i.e.,  at  one  ionization  time  con- 
stant after  ionization  has  been  assumed  to  begin,  we  may  write  the  field-line 
content  as; 


N(r,t) 


(3.6) 
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I(r,  t)  represents  a normalized  field>line  content.  The  normalization  is  chosen 
so  that 

<30 

^yi(r,t)2rdr  = |^1  - exp(-(t  - t.)/T.)J  . (3.8) 

^ 0 

which  approaches  one  as  t -»  <»  . 

We  have  carried  out  the  integration  indicated  in  Eq.  (3.7)  in  order  to 
determine  the  resulting  distribution  of  ions  on  field  lines. 

The  four  thin  solid  curves  in  Figure  3 . 1 show  the  normalized  final 
(as  t -» <»)  field-line  content,  I(r, «),  for  several  different  delay  times,  t.  , 
before  the  onset  of  photoionization.  As  t.  increases,  the  radius  of  the  neutral 
cloud  is  larger  at  the  time  that  photoionization  begins.  As  to  be  expected,  the 
ionization  is  less  peaked  at  the  center  and  is  broader  for  larger  t^ . Note  that 
each  of  the  curves  is  scaled  to  the  radius,  a,  of  the  neutral  cloud  at  t = t.  + 
thus,  the  actual  spread  in  kilometers  for  increasing  t.  is  greater  than  indicated. 

The  heavy  curve  shows  the  assumed  Gaussian  shape  of  the  neutral  cloud 
on  the  same  scale.  This  is  the  sh:q)e  that  the  ion  content  would  have  if  all  of  the 
ionization  were  created  in  a time  small  compared  to  the  expansion  time  of  the 
neutral  cloud,  i.  e. , as  if  all  of  the  ionization  were  deposited  instantaneously  at 
t = • Figure  3. 1 indicates  that  the  scale  size  of  the  ion  cloud  transverse 

to  the  magnetic  field  is  considerably  narrower  than  the  scale  size  of  the  neutral 
cloud.  It  is  also  considerably  narrower  than  that  used  in  the  previous  modeling.  ^ 
The  area  under  each  of  the  five  curves  multiplied  by  2 r/a  is  equal  to  one  indi- 
cating that  each  curve  represents  the  same  total  inventory  of  ions. 
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NORMALIZED  FIELD-LINE  CONTENT,  I(r,oo) 
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Figure  3.1.  Normalized  field-line  ion  content  deposited  by  the  expanding 

neutral  cloud  as  a function  of  radius  transverse  to  the  magnetic 
field  line.  Four  curves  are  shovm  for  photoionization  beginning 
at  four  different  times,  t.  , after  release. 


We  will  make  use  of  this  ion  depositton  model  in  order  to  derive 

parameters  appropriate  for  the  Spruce  ion  cloud  in  Section  5.  We  wiil  mention 

ere  that  the  total  field  line  content  at  the  center  of  the  ion  cioud  according  io 
model  A is  2.4  x ® 


4.  MOTION  AND  DESCENT  OF  THE  ION  CLOUD 


In  this  section  we  discuss  the  motion  of  the  ion  cloud  in  the  vertical 
direction  and  derive  the  distribution  of  ionization  in  the  ion  cloud  along  mag- 
netic field  lines.  We  find  that  the  simple  but  useful  model  described  here 
allows  us  to  obtain  the  rate  at  which  the  ion  cloud  settles  under  the  influence 
of  gravity  and  neutral  winds,  and  the  change  in  sh£^e  of  the  distribution  of 
ionization  parallel  to  the  magnetic  field.  We  find  that  the  ion  cloud  evolves 
into  a characteristic  shape  rather  quickly  and  maintains  this  shape  as  it 
descends  through  the  upper  atmosphere.  The  model  includes  the  response  of 
the  ion  cloud  to  the  force  of  gravity,  spreading  of  ionization  along  the  mag- 
netic field  lines  due  to  ambipolar  diffusion,  and  the  flattening  of  the  bottom 
side  of  the  ion  cloud  as  it  descends  into  the  denser  neutral  atmosphere . 

From  this  model  we  determine  the  peak  electron  density  in  the  ion 
cloud  as  a function  of  time  as  well  as  describe  the  change  in  the  character- 
istic width  of  the  peak  of  the  electron  density.*  Furthermore,  we  shall  calcu- 
late the  change  in  the  Pedersen  conductivity  of  the  ion  cloud  as  it  descends. 
Lastly,  we  shall  compare  the  results  of  the  analytic  solution  obtained  here 
with  the  results  of  the  numerical  solution  obtained  by  Chavin  and  Kilb. 

The  deformations  that  take  place  in  the  ion  cloud  contours  in  the 
plane  perpendicular  to  the  magnetic  field  will  not  be  discussed  in  this  section. 
These  deformations,  ultimately  resulting  in  the  development  of  striations, 
have  been  the  subject  of  other  analytical  and  numerical  studies.  Neither  will 
we  be  concerned  here  with  the  factors  that  lead  to  stabilization  of  short  wave- 
lengths or  mechanisms  that  control  the  diffusion  of  the  plasma  cloud  perpen- 
dicular to  the  magnetic  field. 

*Eilectron  concentration  and  ion  concentration  are  e<{ual  to  each  other  and 
will  be  used  interchangeably  throughout  this  report.  Electron  density  or  ion 
density  is  sometimes  used  to  mean  the  same  quantity. 
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4. 1 BASIC  MODEL  AND  EQUATIONS 

For  our  model  we  adopt  a number  of  simplifying  assumptions.  We 
assume  that  the  ions  and  electrons  are  isothermal  with  temperatures  T.  and 
respectively . We  take  the  magnetic  field,  3 , as  uniform  and  constant 
but  making  a dip  angle,  8 , with  the  horizontal  direction.  We  assume  that  the 
number  density  of  the  neutral  atmosphere,  n^  , varies  exponentially  with 
height,  h , with  a constant  scale  height,  , i.e.,  n^  = n^^exp|  -|h-h^j/Hj^j. 
We  also  assume  that  the  electrons  are  completely  collisionless,  undergoing  no 
collisions  with  either  neutrals  or  ions. 

We  adopt  a coordinate  system  oriented  with  the  magnetic  field  which 
we  take  to  be  pointing  downward  and  northward.  Positive  x will  be  towards 
the  east.  The  z -direction  is  anti -parallel  to  S and  positive  y points  toward 
the  north  and  upwards.  The  acceleration  due  to  gravity,  g , is  vertically  down- 
wards and  has  components  g^  = - g sin  0 and  g^  = - g cos  6 . 

By  setting  S In  Eq.  (1.1)  equal  to  0,  we  write  the  ion  continuity 
equation  as 


Sn 

at 


+ 


-* 

7 • nv. 

X 11 


+ 


0 


(4.1) 


where  n is  the  plasma  density  and  the  ion  velocity  v.  is  given  explicitly  in 
terms  of  components  anti -parallel  to  S and  perpendicular  to  B . The  equa- 
tions of  motion  for  the  ions  and  the  electrons  are 


0 = e(S  + Vj  X S)  - T.7'tn(n)  - i/m(vj  - v^j  + Mg  , (4.2a) 

0 = - e(E  + Vg  X S)  - T^7<n(n)  , (4.2b) 

where  S is  the  electrostatic  field,  v^  is  the  neutral  wind  velocity,  M is  the 
ion  mass,  u is  the  ion-neutral  collision  frequency,  and  the  temperatures  are 
given  in  energy  units.  We  have  neglected  the  inertial  terms  because  the  devel- 
opment of  the  cloud  takes  place  on  a time  scale  typically  measured  in  minutes 
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which  is  long  compared  to  the  collision  time  r = 1/v  which  has  typical  values 
of  0. 1 - 2. 0 s depending  on  altitude.  Our  model  treats  the  electrons  as  being 
completely  collisionless  and  massless. 

It  is  convenient  to  explicitly  separate  the  ambipolar  part  of  the 
electrostatic  field,  (TgV-tn(n|/e  , from  the  remaining  part,  , which 
remains  constant  along  a magnetic  field  line  (as  we  shall  see  below) . Thus, 
we  write  the  electrostatic  field  as 

T 

E = - — 7-tn(n)  . (4.3) 

Note  that  since  v x ^ = 0 , we  also  have  v x E * = 0 . If  we  define 
u ■ v^  + g r , we  may  rewrite  the  ion  and  electron  momentum  equations  as 

0 = e(E'  + V.  X b)  - 2Tvto(n)  - vM^v.  - u)  , (4.4a) 

0 = -e(E*  + Vg  X S)  , (4.4b) 

where  we  have  defined  T as  equal  to  the  average  of  the  ion  and  electron 
temperatures,  T = (t.  + /2  . In  this  form  we  can  see  that  the  effect 

of  gravity  on  the  ion  cloud  is  indestinguishable  from  the  effect  of  a neutral 
wind  vith  velocity  gr  Although  Perkins  et  al.^^  state  that  gravity  acts 
like  an  electric  field  in  the  potential  equation,  it  actually  appears  precisely  as 
a neutral  wind  because  it  creates  a body  force  on  the  plasma  in  the  same  man- 
ner that  collisions  with  neutrals  do  and  exerts  no  force  on  the  electrons  as 
an  electric  field  would. 

The  force  equation  on  the  plasma  cloud  is  obtained  by  adding  Eqs. 
(4.4a  & b)  resulting  in 

7 X 5 = vp  + i/nM^Vj  - u)  (4.5) 
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where  1 = e (vj  - is  the  current  density  and  p = n (t^  + j is  the 
total  plasma  pressure.  In  this  form  we  see  that  the  current  flowing  perpendi- 
cular to  the  magnetic  field  consists  of  two  parts:  a diagmagnetic  part  that  pro- 
vides pressure  balance  transverse  to  the  magnetic  field  and  the  remaining  part 
that  overcomes  the  drag  between  the  ions  and  the  neutrals  (and  the  force  of 
gravity) . It  is  the  latter  part  that  determines  the  slip  velocity  transverse  to 
the  magnetic  field  that  produces  the  deformations  of  the  ion  cloud  leading  to 
backside  steepening  and  striation  formation,  subjects  not  discussed  here. 

The  component  of  Eq.  (4.5)  parallel  to  the  magnetic  field  gives  the 
parallel  ion  velocity 


V. 

iz 


u 

z 


i/M  S z ' ' 


(4.6) 


The  component  of  Eq.  (4.4b)  parallel  to  the  magnetic  field  shows  that  = 0 . 
This  fact  together  with  the  fact  that  2'  is  curl -free,  v x 2'  = 0,  shows 
that  S * has  components  only  perpendicular  to  the  magnetic  field  and  these 
components  are  functions  only  of  the  coordinates  transverse  to  the  magnetic 
field  independent  of  z . (Remember  that  z does  not  measure  distance  in  the 
vertical  direction  but  distance  along  the  magnetic  field  line.) 

In  order  to  obtain  the  ion  velocity  perpendicular  to  the  magnetic 
field  from  Eq.  (4.4a)  it  is  convenient  to  introduce  the  vector  F which  has 
only  components  perpendicular  to  the  magnetic  field 


^ , ir  X £ - II  v^to(n) 


(4.7) 


where  b = B/B^=  ~ the  unit  vector  in  the  direction  of  the  magnetic 

field.  By  defining  ( as  the  ratio  of  the  ion -neutral  collision  frequency  to  the 
ion  gyrofrequency,  (I, 


VIA  _ 1 

^ = 71"  07  ’ 


(4.8) 
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we  may  represent  the  solution  of  Eq.  {4.4a)  for  the  perpendicular  component 
of  the  ion  velocity  as 


''li  = — 2-  (J'  b + f f) 

= (?  - t ? X b)  , (4.9) 

1 + c 

where  = g — - ^ 7 to(n)  x b 

B 

= 7-tn(n)x  b . (4.10) 


At  high  altitudes,  above  140  km,  the  parameter  c < < 1 . To  lowest 
order  in  e the  ion  drift  velocity  is  given  by  v . By  recognizing  that 

A xo 

V • vn  X b = 0 and  v x E^  = 0 , the  continuity  equation  for  the  ions,  Eq. 
(4. 1),  can  be  rewritten  as 


in  E'xB 


Vj^  n + 


2T  ?n 
vM  iz 


0 . 


(4.11) 


!We  thus  see  that  in  lowest  order  the  perpendicular  motion  of  the  ion  cloud 

corresponds  to  incompressable  flow  transverse  to  the  magnetic  field. 

The  equation  that  determines  ? and  thus  the  electrostatic  field  2 * 
from  Eq.  (4.7)  is  obtained  by  first  setting  v • 7 = ev  • n^v.  ■ = 0 

which  can  be  written 


1+  f 


n ? = 


1+  c' 
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n F X b — -r*“ 
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a z ^2 


(4.12) 
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The  first  term  on  the  right -hand  >side  represents  the  divergence  of  the  Hall 
contribution  to  the  current  which  is  less  than  the  Pedersen  current  given  on  the 
left -hand -side.  An  ^proximate  equation  for  which  depends  only  on  coor- 
dinates transverse  to  the  magnetic  field  is  obtained  by  integrating  Eq.  (4.12) 
along  the  magnetic  field  line  and  neglecting  the  Hall  contribution  to  the  current 
which  is  higher  order  in  e . Even  if  the  region  of  integration  extends  to  low 
altitude  where  f ^ 1 , this  approximation  is  still  valid  if  the  field -line -integrated 
Hall  current  is  much  less  than  the  Pedersen  current.  Although  this  approxima- 
tion not  valid  during  the  daytime  or  in  the  auroral  zone  during  auroral  pre- 
cipitation, it  is  valid  after  sunset  at  midlatitudes.  The  pressure  contribution 
to  F is  relatively  small  so  for  simplicity  we  neglect  this  term  for  the  moment. 
This  neglect  has  no  essential  impact  on  the  arguments  which  follow. 

By  carrying  out  the  integration  of  Eq.  (4. 12)  parallel  to  B with 
the  above  approximations,  and  making  the  usual  assumption  that  no  current 
flows  out  of  the  ionosphere  along  magnetic  field  lines,  the  equation  for  2'  can 
be  simply  written 


• r (x,  y)  [E ' + U X = 0 . (4. 13) 

E is  the  field -line -integrated  Pedersen  conductivity  and  IJ  is  a field-line - 
weighted  average  of  the  velocity  field,  u , defined  by 


(4.14) 


(4.15) 


Note  that  although  £ and  are  explicit  functions  of  the  transverse  coordinates, 
they  are  also  functionals  of  the  distribution  of  ionization  along  the  magnetic 
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field  and  thus  of  the  z -coordinate  of  the  center  of  the  cloud,  z . Note  that 
-*  -*  c 

if  u = V + g r where  v = constant  and  if  we  define 
i ni  ni 


v(x, 


then 


-♦ 


U 


g cos  9 
u 


The  set  of  equations  (4. 11),  (4. 13),  (4. 14)  and  (4. 15)  is  the 

fundamental  set  that  controls  the  motion  of  plasma  clouds  in  the  ionosphere. 

11  12  12  14 

This  set  is  equivalent  to  the  standard  equations  used  by  other  workers . ’ ’ ’ 

The  principle  approximations  made  are:  (1)  the  neglect  of  electron  collisions 
whose  principal  influence  is  to  cause  E'  to  have  a slight  variation  along  the 
magnetic  field;  and  (2)  the  neglect  of  Hall  current  and  other  terms  of  higher 
order  in  c . Neither  of  these  effects  has  an  influence  on  either  the  motion  of 
ionization  along  the  magnetic  field  or  on  the  zeroth  order  motion  transverse 
to  the  field  that  are  the  subjects  of  discussion  in  the  following  subsections. 

The  importance  of  these  neglected  effects  on  striation  behavior  and  the  amount 
of  sophistication  of  treatment  needed  to  incorporate  them  is  an  unsettled  ques- 
tion and  an  area  of  active  investigation  beyond  the  scope  of  this  report. 

These  equations  do  incorporate  additional  features  not  ordinarily 
included  that  do  have  a bearing  on  the  subject  of  the  behavior  and  motion  of 
the  ion  cloud.  They  include  a neutral  atmosphere  whose  density  is  a function 
of  altitude  and  a neutral  wind  velocity  field  that  is  not  spatially  uniform.  The 
effect  of  gravity  is  included  and  has  been  shown  to  have  the  same  effect  as  a 
downward  neutral  wind  whose  magnitude  is  inversely  proportional  to  the  neu- 
tral density . In  the  next  subsection  we  discuss  the  motion  of  the  ion  cloud 
that  results  from  these  equations. 
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4.2 


MOTION  OF  THE  ION  CLOUD 


The  qualitative  properties  of  the  solution  of  Eq.  (4. 13)  are  well 
known.  For  small  plasma  clouds  that  make  a very  small  chaise  in  the 
Pedersen  conductivity,  Z , the  quantity  in  brackets  is  essentially  constant 
and  E has  essentially  the  value  of  the  ambient  electric  field  Eo  • 0*^  the 
other  hand,  high-density  plasma  clouds  that  cause  significant  change  in  the 
field-line -integrated  Pedersen  conductivity  causes  the  term  in  brackets  to  be 
nearly  zero  for  those  magnetic  field  lines  that  pass  through  the  center  of  the 
ion  cloud.  Thus  E'  for  those  field -lines  has  the  approximate  value  - U x B. 

We  may  derive  useful  approximate  e3q>ressions  for  the  velocity  of  the 
center  of  the  ion  cloud  transverse  to  the  m^^etic  field  line  by  separately 
identifying  the  ambient  ionospheric  density  n^  and  the  ion  cloud  density  n^  . 
For  the  field  line  that  passes  through  the  center  of  the  cloud  at  the  coordinates 
X = x^  and  y = y^  , we  write 

^(*c'  Vc)  = . (4.16) 

where  I is  the  field-line -integrated  Pedersen  conductivity  of  the  cloud,  and 
Z the  field -line -integrated  Pedersen  conductivity  of  the  ionosphere.  The 
velocity  of  thecenter  of  the  ion  cloud  transverse  to  the  magnetic  field,  , 
can  be  written 


S'  K, 


'c’ 


B' 
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(4-17) 


where  tj  is  the  coupling  coefficient  of  the  ion  cloud  to  the  field -line-weighted 
neutral  wind  U.  i)  is  a function  of  Z*^/Z^  that  has  the  following  properties: 
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Z^  < < Z^ 
Z*^  > > Z^ 


(4.18) 


Note  that  this  definition  of  a coupling  parameter  differs  from  previous 
descriptions  in  terms  of  the  parameter  f = 1 - 17  . 


I 


If  we  further  define 


Ua  = 


/.T7“- 

/ 

J 1 + 


dz 


a = 


“■■.c  ' (’■c-  >'c-  ^c)  • 

then  we  may  define  U as  having  the  approximate  value 

^ ^c)  - “ ®a  * rfsPic  - "a) 


(4.19) 


For  a cylindrically  symmetric  model  ion  cloud  with  a constant  den- 
sity inside  a finite  radius  and  zero  density  outside,  T7(ot)  * has  the  following 
14 

functional  form, 


f7(«)  s a/(2  + a)  . (4.20) 

15  16 

other  studies  of  the  influence  of  shape  on  the  motion  of  ion  clouds  ’ have 
shown  that  this  approximation  to  ri  (a)  slightly  over-estimates  the  actual 
amount  of  coupling  to  the  neutral  wind  for  ion  clouds  with  more  realistic  trans- 
verse density  profiles.  However,  in  terms  of  ri{a)  given  by  Eq.  (4.20),  we 
may  write  an  explicit  expression  for  the  velocity  of  an  ion  cloud  transverse  to 
the  magnetic  field; 
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(1  + 0^(2  + 0$  2 + a g2 


(4.21) 


This  expression  should  be  approximately  valid  throughout  the  range  0 < et  < oo  • 

In  the  appropriate  limits  it  says  that  lovrdensity  clouds  (a  < < 1)  move  per- 
pendicular  to  B with  the  ambient  E„  x B/te  velocity  while  high-density 

A 

clouds  (a  > > 1)  are  carried  along  with  the  neutral  wind  at  the  location  of  the 
cloud. 

The  motion  parallel  to  the  magnetic  field  is  precisely 

u = u (x  , y ; z \ . Hence  the  motion  of  the  cloud  is  obtained  by  combin- 
zc  zy  c c c)  i 

ing  the  perpendicular  and  parallel  motions.  For  instance,  we  will  derive 
expressions  for  the  downward  motion  of  an  ion  cloud  in  the  earth  frame  in 
terms  of  the  northward,  v^^  , and  downward,  Vj^  , components  of  the  neutral 
wind.  Eastward  components  of  the  neutral  wind  do  not  affect  vertical  motion 
of  ion  clouds.  Since 

“j.  = jvjj  sin  0 - (Vjj  + g t)  cos  ej  ly  (4.22) 

and  denoting  the  y-component  of  U|x^,  y^;  z^j  , (see  Eq.  (4.15))  by  , and 
by  using  Eq.  (4.17),  we  may  write  the  downward  component  of  velocity  of  an 
ion  cloud  as 


= - tj  cos  0 + |^(vjj  + g rj  sin  0 + v^j  cos  0j  sin  0 + (1-tj) 


cos  0 


'j 


where  E is  the  westward  component  of  the  ambient  electric  field, 
w 

For  low-density  clouds,  the  downward  velocity  is  obtained  by  setting 
T)  = 0 in  Eq.  (4.23)  which  can  be  written 


cD 


~ u sin  0 + 
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cos  0 


(4.24) 
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i.e.,  the  cloud  motion  is  simply  the  cdmponent  of  u parallel  to  B plus  the 
E X B motion  transverse  to  B . A low-density  cloud  slides  up  and  down  the 
field  lines  in  response  to  the  neutral  forces  parallel  to  B and  crosses  field 
lines  according  to  E x S motion.  Thus  northward  and  downward  neutral 
winds,  gravity,  and  westward  electric  fields  all  contribute  to  downward 
motion  of  low -density  clouds. 

For  high-density  clouds,  ri  = 1 and  U is  the  y -component  of  u 

W X 

given  by  Eq.  (4.22).  The  expression  (4.23)  reduces  to 

'^cD  ~ ^ ST  (4.25) 

i.e.,  the  ion  cloud  descends  as  thoi^h  there  were  no  electric  and  magnetic 
fields  present.  One  word  of  caution  must  be  raised  at  this  point.  Although 
high -density  ion  clouds  move  as  though  there  were  no  magnetic  and  electric 
field  present,  it  must  not  be  assumed  that  the  ionization  behaves  as  though  it 
were  unchained.  If  Eq.  (4.11)  is  integrated  parallel  to  the  magnetic  field,  the 
resulting  equation  shows  that  the  total  field-line -integrated  ionization  content 
is  a conserved  quantity . Every  piece  of  ionization  on  a given  field  line  moves 
across  magnetic  field  lines  at  the  same  velocity  even  in  the  limit  of  high- 
density  clouds. 

4. 3 VERTICAL  DISTRIBUTION  OF  IONIZATION 

We  now  turn  our  attention  to  determining  the  vertical  motion  of  the 
ion  cloud  and  its  change  of  shape  along  the  magnetic  field  line.  We  shall 
restrict  ourselves  to  the  special  but  reasonable  case  that  the  neutral  wind 
field  depends  only  on  altitude 

h - h^  = y cos  6 + z sin  9 (4.26) 

where  h^  is  the  initial  altitude  of  the  cloud  and  the  origin  of  the  x,  y,  z 
coordinate  system  is  at  the  initial  location  of  the  center  of  the  cloud.  The 
altitude  of  the  center  of  the  cloud,  h (t)  , can  be  calculated  as  a function  of 
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time, 


t = 


dh 

J v.dII') 


(4.27) 


where  is  the  altitude -dependent  downward  velocity.  The  velocity  of  the 

ion  cloud  in  the  y -direction,  V (h  (t)).  as  given  by  combining  Eqs.  (4.22) 

cy  c 

and  (4.19)  with  (4.17),  is  now  an  explicit  function  of  t which  can  be  integrated 
to  obtain  the  y -coordinate  of  the  field  line  passing  through  the  center  of  the  ion 
cloud, 

V*' 
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We  will  now  show  that  the  changes  in  density  that  occur  along 
nonvertical  magnetic  field  lines  for  a particular  class  of  altitude -dependent 
neutral  winds  can  be  described  by  an  equation  valid  for  vertical  magnetic  fields 
and  no  neutral  wind.  The  distribution  of  ionization  along  the  central  magnetic 
field  line  is  given  by  the  function 

n^(z,  t)  = n(x^,  y^,  z,  t)  , (4.29) 

I where  and  y^  are  the  coordinates  of  the  field  line  passing  through  the 

! center  of  the  cloud.  The  equation  for  n (z,  t)  is  obtained  by  transforming 

|j  ^ -t 

Eq.  (4. 11)  to  the  frame  of  reference  that  is  moving  with  the  cloud  velocity 
I given  by  Eq.  (4. 17)  transverse  to  the  magnetic  field.  This  transformation 

I allows  us  to  write  the  ion  continuity  equation  in  the  form 
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(4.30) 


where  u and  v are  functions  of  the  altitude  along  the  field  line  passing 

ZC  c 

through  the  center  of  the  cloud. 


50 


(4.31) 


h = (t)  cos  6 + z sin  0 . 

The  ion -neutral  collision  frequency  is  proportional  to  the  atmospheric 
density  which  we  will  assume  varies  exponentially  with  a constant  atmospheric 
scale -height,  , 

u = y^exp(-{h-ho)/H_,)  . (4.32) 

Furthermore,  we  assume  that  the  component  of  u along  the  magnetic  field 
has  the  explicit  altitude  dependence 

“z  = - - «z2  Mi”  - "oK)  • 

is  a constant  downward  drift  parallel  to  the  magnetic  field  while  the  sec- 
ond term  corresponds  to  a constant  downward  flux  of  neutral  particles.  We  also 
take  the  transverse  component  of  the  ion  cloud  velocity  to  have  the  same  form; 

''ey  = -V  ■ 

With  the  definitions  = Vyj  cos  0 + sin  0 and  Vjjg  “ ^y2  ® 

+ U^2  ® downward  velocity  of  the  center  of  the  cloud  has  the  same 

form, 

''oD  = ''dI  ♦ ''d2'=^(('’-''oK) 

We  may  now  explicitly  carry  out  the  integrations  given  in  Eqs.  (4.27) 
and  (4.28)  which  yield 
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(4.37) 


Note  that  for  times  longer  than  it  takes  to  drift  a scale-height,  , the 

cloud  descends  at  constant  velocity,  and  has  covered  a distance  only 

H^-tn  (1  + greater  than  the  distance  . If  there  is  no 

vertical  constant  drift,  Vj^j  = 0 , then  Eq.  (4.36)  reduces  to 

•■c®  -“o  . (4.38) 


Note  that  in  this  case,  the  downward  velocity 

dhc(t) 
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dt 
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n 


(4.39) 


at  first  is  as  to  be  expected  but  decreases  only  as  t"^  and  not 

exponentially  at  late  time. 

We  now  introduce  a new  coordinate  z*  = z + and  transform 

Eq.  (4.30)  to  a frame  drifting  down  the  magnetic  field  line  with  speed  . 
Set  n^(z,  t)  = ff(z*,t),  and  obtain 
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where  = 2T/v^M  and 


d*(t)  = 


U^jtsinS  - y^(t)  cos  9 


H_ 


._L  ^D2 


(4.41) 


d’"  (t)  is  the  vertical  distance  downward  measured  in  scale -height  units 
that  the  cloud  has  moved  due  to  the  constant  drift  along  field  lines 
plus  the  transverse  motion  across  field  lines.  Equation  (4.40)  describes  the 
evolution  of  the  ion  cloud.  If  = 0 0 = 0 corresponding  to 

horizontal  magnetic  field  lines,  it  describes  diffusion  of  the  cloud  with  a time- 
dependent  diffusion  coefficient  as  the  cloud  drifts  downward  across  field  lines. 
If  U^2  ^ ® is  siso  a time -dependent  drift  along  field  lines  which  does 

not  influence  the  change  in  shape.  If  sin  6 / 0 , the  effective  diffusion  coef- 
ficient is  also  a function  of  position  along  the  field -line. 

Equation  (4.40)  can  be  brought  into  the  simple  form 


an*  d 


r z*/ai 


+ an*l 


where  z*  = z'  sin  e/H^  , n(z',  t)  = n*(z*,  t*),  a 


(4.42) 


sin  0 and 

Ztt  n o 


2 ^ 
D Sin  0 r 


n o 


^z2  ®^D2 


- 1 . (4.43) 


n*(z*,t*)  gives  the  density  distribution  as  a function  of  altitude  measured  in 

2 2 

scale  height  units  . As  t -*  0 , t*  -♦  D^t  sin  ®/H^  and  as  t -*  <» , 
t*  •*  constant  unless  = 0 . If  gravity  alone  acts  to  pull  particles  down 

field  lines,  then  = S sin  B/i/^  and  we  can  write 


Tn  M 


(4.44) 


where  is  the  barium  plasma  scale  height,  and  are  the  atmospheric 
temperature  and  the  mean  molecular  weight.  For  electrons  only  a little  hotter 
than  ions  and  molecular  weights  around  20,  a has  a value  around  3. 

Equation  (4.42)  describes  the  descent  of  neutral  particles  in  an 

17 

exponential  atmosphere  and  its  solution  was  given  by  Yu  and  Klein  in  the  form 

nnz*,t*)  = .C  ^a+l^^|’_  (4  45) 

where  (,  = exp  (- z*/2)  and  I . is  the  modified  Bessel  function  of  order 
a-1  . The  normalization  constant  is  determined  by  integrating  the  cloud  den- 
sity along  the  field  line  and  designating  the  maximum  total  field  line  content 

of  the  cloud  by  N , 
c 

00  H ^ H 

Nc  = = -SiTT «•«) 


N sin  e 


so  that  C = 


H. 


n 


At  early  and  late  times,  the  approximations  to  the  Bessel  function 

18 

valid  for  large  and  small  values,  respectively,  of  its  argument  can  be  ^plied 
to  give  the  behavior  of  the  density  at  early  and  late  times: 


n(z*,  t*) 


(4.47) 


(4.48) 


where  for  the  sceond  expression  t*  » 1 and  the  maximum  value  occurs  at 
C = = at*  . At  early  time,  the  solution  given  by  Eq.  (4.47)  corresponds 

to  diffusion  of  a point  release  in  a constant  density  atmosphere  and  hence  is 
applicable  to  barium  releases  from  about  a minute  after  release  after  the  ion 
deposition  phase  discussed  in  the  preceeding  section  has  finished. 


The  asymptotic  form  valid  at  late  time  given  by  Eq.  (4. 48)  can  be 
compared  to  another  exact  solution  of  Eq.  (4. 42),  n , which  is  valid  for  all 
values  of  t*  but  satisfies  an  initial  condition  different  from  that  satisfied  by 
Eq.  (4.45); 


where 


m e'^m  = 1 + at* 
m 


"m))-  ‘]) 

(4.49) 

(4.50) 


55 


2 — 

is  the  value  of  4 ^t  which  n reaches  its  maximum  value.  The  only 

18 

£^proximation  in  the  last  expression  is  in  using  Stirling's  formula  for  r(a)  . 
By  using  Eq.  (4.43)  for  t*  and  recognizing  that 

d^(t)  - zj;^(t)  = (h^  - h^(t))/H^  , (4.51) 


we  can  write  the  solution  n as  an  explicit  function  of  altitude 


n(z*,  t) 


The  solution  n corresponds  to  a density  distribution  of  constant  shape  with  its 
maximum  at  h = h^(t)  given  by  Eq.  (4.36).  Near  the  peak,  the  density  has  a 
Gaussian  shape, 


5(h-.h^,t) 


N sin  e (h  - h (t))* 

exp k 


(4.53) 


but  decays  exponentially,  n exp  ^-|h  - h^j  altitude  with  a 

scale-height  • It  decays  more  rapidly  than  a Gaussian  fall -off  on  the  lower 

side. 


The  density  distribution  given  n in  Eq.  (4.49)  is  identical  to  that 

given  by  the  asymptotic  distribution  in  Eq.  (4.48)  to  order  1/t*  . The  location 

2 

of  the  peak  given  by  = 1 + at’*'  is  accurate  but  not  exact  and  has  the  correct 
values  for  both  small  and  large  t’*'  . Numerical  evaluation  shows  that  the  shape  of 
n(z*,  t)  is  well  approximated  by  n(z*,  t)  for  at’*'  ?>  2 . 

“ 19 

The  shape  n is  the  same  as  that  discussed  by  Rao  in  connection  with 

barium  clouds  for  the  case  11^2  = g r sin  6 and  no  neutral  winds.  We  have 

derived  the  fully  time -dependent  solution  for  the  density  distribution  including 
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altitude -dependent  winds  and  magnetic  field  of  arbitrary  dip.  Before  giving 
examples  for  specific  cases  we  make  some  general  observations . 

The  factors  that  influence  the  shape  of  the  density  distribution  are 
the  values  of  t*  and 


^z2  = '\2  ^ ' 

that  is,  gravity  and  an  altitude-dependent  component  of  the  neutral  wind 
parallel  to  the  magnetic  field.  The  value  of  t*  is  strongly  influenced  by  the 
downward  displacement,  d*  , of  the  cloud.  If  d*  is  small,  an  under- 
estimate of  the  time,  t , it  takes  the  cloud  to  evolve  to  the  asymptotic 

A 

shape  given  by  n is  obtained  by  setting  the  early -time  Gaussian  width  (see 

2 2 i 

Eq.  (4.47) ) equal  to  the  asymptotic  value  (2  H /a  sin  6)  . This  pro- 

2 2 

cedure  results  in  t H /2  a D sin  6 which  is  equivalent  to  (see  Eq. 
a n o 

(4.43) ) a t*  ~ 0.5  or  about  1/4  of  the  time  indicated  by  numerical  eval- 
uation. After  4 t , the  electron  density  in  the  vicinity  of  the  peak  is  pro- 
“ 1 

portional  to  (U^g/s  ^ Gaussian  scale  length  along  the  mag- 

netic field  is  H*  = (2  H H g t/U  - sin  0)^  . 

n c 

However,  downward  motion  of  the  cloud  transverse  to  the  magnetic 
field  and  a constant  drift  parallel  to  the  magnetic  field  each  have  a strong 
influence  on  the  time  development  of  the  cloud.  If  = 0 , no  constant 
downward  drift  of  the  cloud  but  only  altitude^ependent  downward  drift  across 
field  lines,  then  it  takes  longer  than  4t  to  reach  the  asymptotic  shape  n . 

A 

If  there  is  a constant  downward  drift,  0 , then  inspection  of  Eq.  (4.43) 

shows  that  t*'  approaches  a constant  value,  t|  . Generally  at|  ~ 1 but  if 
>>  V^2  ’ value  of  at^  can  be  much  less  than  1 . In  this  case  of 
rapid  descent  of  the  ion  cloud,  its  final  shape  would  remain  Gaussian -like 
with  a narrower,  higher -valued  peak  than  that  predicted  by  n . 
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4.4 


DENSITY  PROFILES  FOR  SPEQFIC  CASES 


The  effect  of  the  neutral  wind  on  the  time  development  of  an  ion  cloud 
is  discussed  here  for  several  specific  cases.  In  particular,  we  show  the  effects 
due  to  altitude -dependent  and  altitude -independent  neutral  winds,  downward  and 
upward  neutral  winds,  and  the  response  of  both  small  and  large  barium  releases. 

The  formulas  derived  in  the'  last  subsection  depend  on  an  atmospheric 
density  model  in  which  the  ion -neutral  collision  frequency  varies  exponentially 
with  altitude  according  to  Eq.  (4.32).  An  expression  for  the  collision  time, 

T = 1/v  , can  be  obtained  from  measurements  of  the  mobility  of  barium  ions  in 
nitrogen  gai^®  resulting  in 


T 


8.6  X 10 


15 


n 


(4.54) 


where  n^^  is  an  effective  neutral  concentration  in  particles/m 


3 


n 


n 


+ n^  + 0.8n. 


(4.55) 


In  deriving  Eq.  (4.54)  it  has  been  assumed  that  the  collision  cross  section 
varies  inversely  as  the  square  root  of  the  temperature  and  that  the  collision 
cross  section  for  oxygen  molecules  and  atoms  is  the  same  as  it  is  for  nitrogen 
molecules.  The  fact  that  the  lighter  oxygen  atoms  are  less  efficient  in  stopping 
a heavier  barium  ion  is  expressed  by  the  appropriate  Langevin  factor, 

(‘  " (>  * ' »'8  ■ 


where  and  Mq  are  the  masses  of  the  nitrogen  molecule  and  oxygen  atom 

2 

respectively . 


In  order  that  a comparison  of  the  results  of  our  analytic  calculations 
can  be  made  with  the  numerical  values  obtained  by  Chavin  and  Kilb^^when  they 
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did  a simulation  of  event  Anne,  we  have  used  their  atmospheric  density  model. 
The  solid  curve  in  Figure  4.1  shows  the  effective  density,  n,  obtained  using 
their  values.  The  dashed  line  shows  an  ej^nential  fit  to  the  height  variation 
of  the  effective  density  resulting  in  an  effective  scale  height  of  22  km.  This 
scale  height  is  somewhat  less  than  the  atmospheric  density  scale  height  because 
the  concentration  of  oxygen  atoms,  vdiich  are  dominant  above  150  kilometers 
altitude,  increases  less  rapidly  than  the  mass  density  at  lower  altitudes  as  a 
result  of  chemical  reactions.  By  combining  the  above  numbers  we  obtain  the 
desired  expression  for  the  altitude -dependence  of  the  product  g r , 

gT  = gT^  exp  ( (h  - 180) /22)  (4.56) 

where  = 7.7  m/s  . 

General  expressions  for  the  ion  cloud  velocity  perpendicular  to  the 

magnetic  field  and  downward  are  given  by  Eqs.  (4.21)  and  (4.23),  respectively. 

If  u and  have  an  altitude -dependence  in  the  forms  given  by  Eqs.  (4.33) 

Ci 

and  (4.34),  we  find  that  the  general  form  for  the  downward  motion  of  the  ion 
cloud,  Eq.  (4.23),  also  has  the  form  assumed  in  Eq.  (4.35).  For  the  specific 
cases  of  low-density  and  high-density  clouds,  the  downward  velocities  reduce 
to  the  simple  e7q}ressions  given  in  Eqs.  (4.24)  and  (4.25),  respectively  and 
are  expressable  in  the  form  given  by  Eq.  (4.35).  Note  particularly  that  for 
low -density  clouds,  the  presence  of  a westward  component  of  the  ambient 
electric  field  has  the  same  effect  as  a constant  downward  neutral  wind.  The 
westward  directed  electric  field  gives  rise  predominantly  to  a southward 
electric  field  drift  as  well  as  to  a small  downward  component. 

Table  4. 1 provides  a listing  of  the  parameters  used  in  order  to 
illustrate  the  effects  of  different  vertical  neutral  winds.  All  clouds  start  out 
initially  with  a Gaussian  shape  drifting  with  an  ^propriate  constant  velocity. 

It  is  only  at  late  times  that  the  effect  of  different  neutral  winds  becomes  appar- 
ent. Case  1 corresponds  to  no  neutral  wind  with  the  ion  cloud  just  settling 
under  the  influence  of  gravity  alone.  Cases  2,  3 and  4 provide  examples  of 


Table  4.1.  Cases  illustrated. 
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The  symbol  C represents  a constant,  uniform  drift. 
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Figure  4.1.  The  effective  number  density  v&  altitude  derived  from  the 
atmospheric  density  model  used  by  Chavin  and  Kilb.  The 
solid  curve  shows  the  calculated  values  and  the  dashed 
curve  shows  the  exponential  fit  to  the  model  that  is  a good 
approximation  between  150  and  185  kilometers  altitude. 


different  downward  neutral  winds.  In  Case  2 the  velocity  is  assumed  to  vary 
exponentially  with  the  form  given  by  Eq.  (4.35).  Cases  3 and  4 correspond  to 
moderate  and  large  downward  constant  neutral  winds.  The  resulting  values  of 
^Dl’  ^D2  ^ (defined  after  Eq.  (4.42))  for  these  cases  are  given  in  the 

fifth,  sixth  and  seventh  columns,  respectively,  of  Table  4.1. 

Cases  5 and  6 correspond  to  upward  neutral  winds  increasing 
exponentially  in  the  former  case  and  remaining  constant  in  the  latter  but  with 
the  same  magnitudes  as  in  Cases  2 and  3.  For  comparison.  Cases  7,  8 and  9 
correspond  to  the  motion  of  small  ion  clouds  where  the  nonvertical  magnetic 
field  inhibits  the  downward  motion  of  the  cloud.  The  specific  values  are  chosen 
to  correspond  to  Cases  1,  2 and  3 respectively.  We  note  that  the  electric  field 
does  not  influence  the  motion  of  large  ion  clouds  which  move  with  the  neutral 
wind.  An  electric  field  in  a nonvertical  magnetic  field,  however,  has  an  equiv- 
alent effect  on  a small  ion  cloud  that  a constant  neutral  wind  has  on  a large  ion 
cloud. 

The  expression  (4.36)  for  the  altitude  of  the  center  of  the  ion  cloud  has 
been  evaluated  for  each  of  the  9 cases.  The  results  are  shown  in  Figure  4.2.  In 
Case  2 the  cloud  behaves  as  though  there  were  an  effective  gravity  1.67  g.  The 
cloud  reaches  the  same  altitude  as  in  case  1 in  three -fifths  of  the  time.  Cases 
3 and  4 show  the  obvious  result  that  a large  downward  constant  wind  is  more 
effective  in  bringing  the  cloud  to  lower  altitudes.  Case  5 corresponds  to  an 
exponentially  varying  upward  neutral  wind  producing  motion  equivalent  to  a 
value  of  gravity  equal  to  0.33  g.  It  reaches  the  same  altitude  as  the  cloud  in 
Case  1 in  a time  three  times  as  long  In  the  case  of  a constant  upward  neutral 
wind  as  in  Case  6,  the  peak  density  In  the  cloud  moves  to  an  altitude  at  which 
g r is  equal  to  the  upward  velocity.  Cases  7,  8 and  9 show  that  a small  ion 
cloud  does  not  descend  quite  as  rapidly  as  the  corresponding  large  clouds 
shown  in  Cases  1,  2 and  3.  In  general,  the  peak  density  in  the  ion  cloud  will 
not  rise  to  an  altitude  higher  than  the  release  altitude  unless  the  upward  neutral 
wind  has  a value  greater  than  g at  the  altitude  of  release. 
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The  altitudes  of  the  peak  electron  density  obtained  by  Chavin  and  Kilb 
in  their  numerical  simulation  of  Event  Anne  are  shown  plotted  by  X in  Figure 
4.2.  In  their  calculations  they  made  the  equivalent  assumptions  of  a large  ion 
cloud  and  no  electric  fields  or  neutral  winds.  The  agreement  with  the  Case  1 
curve  shows  that  Eq.  (4.36)  is  an  appropriate  description  of  the  time -dependence 
of  the  altitude.  Subsequent  computer  runs  by  Chavin  and  Kilb  in  which  horizontal 
neutral  winds  were  included  showed  that  they  had  no  effect  on  either  the  altitude 
or  density  profile  of  the  ion  cloud  which  is  consistent  with  the  results  for  large 
ion  clouds  as  discussed  in  subsection  4.2. 

We  now  turn  our  attention  to  a discussion  of  the  shape  of  the  electron 
density  distribution  as  a function  of  altitude.  The  next  to  last  column  of  Table 
4.1  lists  the  lapse  time  since  release  for  the  cases  illustrated  in  the  figures 
listed  in  the  last  column.  Figure  4.3  shows  the  time  development  of  a large 
barium  cloud  under  the  influence  of  gravity  alone,  Case  1.  The  settling  due  to 
gravity  and  the  approach  to  *he  asymptotic  shape  given  by  Eq.  (4.52)  is  ^parent. 
The  density  distribution  at  1.5  hours  after  release,  Case  la,  is  shown  with  a 
bold  line  as  this  curve  will  be  repeated  in  the  subsequent  figures  when  we  com- 
pare it  to  conditions  of  various  neutral  winds.  The  time  lapse  of  1.5  hours 
since  release  corresponds  to  t ~ 3 t and  the  bold  curve  is  only  slightly  nar- 
rower  with  a 5%  larger  peak  density  than  the  asymptotic  shape  shown  at  t = 3 
hours. 

We  have  compared  the  shapes  of  the  density  profiles  obtained  from  the 

analytic  function  given  in  Eq.  (4.45)  and  shown  in  Figure  4.3  with  the  results  of 

the  density  profiles  obtained  by  Chavin  and  Kilb^^at  corresponding  tidies.  Because 

their  initial  cloud  had  a far  higher  electron  content  than  our  cloud  which  has  a 
17  -2 

content  of  2.4  x 10  m , we  have  arbitrarily  rescaled  their  profiles  to  have 
comparable  ion  content.  In  general  their  profiles  are  broader  than  the  analytic 
shapes  as  they  employed  a more  realistic  atmospheric  model  with  a value  of 
that  was  20%  larger  than  the  value  we  employed  for  illustrative 

purposes . 
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Figure  4.3.  Electron  density  profiles  as  a function  of  time  in  Case  1 of 
no  neutral  wind.  The  cloud  descends  under  the  influence  of 
gravity  starting  from  an  initial  Gaussian  sh^  at  180  kilo- 
meters and  evolves  to  the  asymptotic  shape  given  by  the 
lowest  curve.  The  data  points  shown  are  from  the  numer- 
ical calculations  of  the  simulation  of  the  Event  Anne  by 
Chavin  and  Kilb. 
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Simulation  results  for  the  two  later  times  are  shown  by  O and  X 
respectively.  Their  larger  value  of  results  in  their  less  rapid  fall -off  of 
density  on  the  high  altitude  side.  It  also  contributes  to  the  lowering  of  the  peak 
electron  density  although  the  precise  fit  in  that  region  is  arbitrary.  The  flatten- 
ing on  the  bottom  side  is  due  in  part  to  the  decreasing  scale  height  in  their  atmos- 
pheric model  (see  Figure  4.1)  and  in  part  results  from  boundary  effects  due  to 
the  approach  to  their  low  altitude  limit.  The  excellent  agreement  between  the 
analytic  curves  and  the  numerical  computation  indicates  that  both  are  giving 
equivalent  descriptions  of  the  physics  that  was  assumed  to  govern  the  descent 
of  the  ion  cloud. 

Figure  4.4  shows  the  effect  of  downward  winds  on  the  density  profile 

of  an  ion  cloud.  The  heavy  curve  is  Case  la  from  the  previous  figure  in  which 

no  wind  has  been  acting  and  the  cloud  has  been  settling  under  gravity  alone. 

Case  2 corresponds  to  an  altitude -dependent  downward  wind  with  a magnitude 

of  0.67  gr  at  the  altitude  of  release.  The  cloud  descends  to  a lower  altitude 
o 

and  its  asymptotic  shape  is  narrower  with  a higher  peak  density . Cases  3 and 
4 correspond  to  constant  downward  winds  of  different  magnitude  shown  at  differ- 
ent times  since  release.  Note  that  in  Case  3 the  ion  cloud  has  arrived  at  a 
lower  altitude  than  the  cloud  in  Case  4 even  though  the  product  of  the  constant 
wind  velocity  times  time  is  only  3/4  as  great  as  in  Case  4.  The  reason  is  that 
in  Case  4 the  ion  cloud  descended  so  rapidly  that  the  settling  due  to  gravity 
j rapidly  became  ineffective. 

i 

■ The  reason  for  including  Case  4 is  to  illustrate  the  effect  that  a larger 

downward  wind  has  on  the  asymptotic  shape.  The  asymptotic  shape  shown  by 
Case  4 (t*  « 0.8  t|)  is  considerably  different  from  that  given  by  the  approxi- 
, mate  shape  n in  Eq.  (4.49).  The  downward  motion  has  brought  the  ion  cloud 

! into  a region  of  dense  atmosphere  such  that  the  diffusion  coefficient  is  too 

small  to  allow  the  cloud  to  diffuse  to  the  equilibrium  shape.  If  the  downward 
neutral  wind  were  to  cease  then  the  ion  cloud  shown  in  Case  4 would  expand 
and  its  peak  density  would  be  lowered  and  its  density  profile  would  assume  the 
shape  that  Case  la  has  but  it  would  take  longer  than  1.5  hours. 
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Figure  4.4.  Density  profiles  for  various  downward  neutral  winds. 
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Figure  4.5  shows  the  effect  of  upward  neutral  winds.  In  Case  5 the 
altitude -dependent  neutral  wind  has  reduced  the  effective  gravity  to  a value  of 
0.33  g.  The  result  is  a broader,  lower-valued  peak  that  descends  slowly  in 
altitude.  Case  5 is  shown  for  two  different  times . The  earlier  time  corres- 
ponds to  Case  la  while  at  the  later  time  (three  times  longer)  the  peak  density 
has  been  lowered  to  the  same  altitude  as  the  peak  density  in  Case  la.  Case  6 
corresponds  to  a constant  upward  neutral  wind.  Since  the  diffusion  coefficient 
at  higher  altitude  is  larger,  this  cloud  has  already  achieved  the  asymptotic 
form  that  Case  la  has  almost  approached.  Note  that  although  these  upward 
winds  were  too  small  to  raise  the  peak  electron  density  above  the  release  alti- 
tude, they  do  transport  significant  ionization  to  higher  altitude. 

Figure  4.6  shows  the  profiles  obtained  by  low-conductivity  ion  clouds 
under  the  influence  of  downward  winds.  Cases  7,  8 and  9 correspond  to  Cases 
la,  2 and  3 for  large  ion  clouds.  In  each  case  the  shapes  of  the  low-density 
clouds  are  comparable  to  the  high  -density  clouds  but  they  have  not  descended 
to  quite  as  low  an  altitude.  This  retardation  in  descent  is  due  to  the  nonvertical 
magnetic  field.  The  differences  in  altitude  would  of  course  be  much  larger  for 
magnetic  fields  with  a dip  angle  smaller  than  60* . 

4.5  DISCUSSION 

The  analysis  carried  out  in  this  section  has  been  motivated  by  several 
issues  that  arose  as  a result  of  the  Stress  Test  Series.  The  primary  questions 
that  were  raised  were:  How  quickly  can  ion  clouds  descend  to  low  altitudes? 

Can  the  peak  electron  density  increase  at  late  time  or  as  the  cloud  descends  in 
altitude?  Can  the  bulk  of  the  ionization  be  confined  to  a narrow  altitude  range? 
In  considering  the  answers  to  these  questions,  it  became  £q)parent  that  the 
neutral  winds  and  E x B drifts  would  have  an  appreciable  effect  on  the  vertical 
transport  of  the  ion  cloud. 

We  have  presented  in  this  section  an  analytical  analysis  of  the  effect 
of  neutral  winds  on  the  vertical  transport  of  ionization.  Analysis  has  included 
the  effect  of  a nonvertical  magnetic  field  and  is  applicable  to  ion  clouds 
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Figure  4.5.  Density  profiles  for  various  upward  neutral  winds. 
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Figure  4.6.  Density  profiles  for  low-conductivity  ion  clouds  all  at  1.5 

hours  after  release.  For  comparison  purposes  the  ion  con- 
tent in  the  low-conductivity  clouds  was  kept  at  2.4  x 10^'  m' 
Cases  7,  8 and  9 are  comparable  to  Cases  la,  2 and  3 shown 
in  Figure  4.4. 
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of  arbitrary  density.  Specific  examples  for  the  limits  of  small  and  large  ion 
clouds  have  been  illustrated.  The  comparison  of  the  analytic  results  obtained 
here  with  the  numerical  results  obtained  by  Chavin  and  Kilb^^  shows  that  these 
different  procedures  produce  similar  results.  The  analytic  results  have 
allowed  us  to  investigate  the  influence  of  different  winds. 

We  have  found  that  for  times  of  order  hours  the  ion  cloud  can  descend 
to  quite  low  altitudes  with  only  moderate  downward  neutral  winds.  We  have 
found  that  the  peak  electron  density  can  remain  high  and  narrow  in  altitude 
only  if  there  is  a rapid  downward  drift.  Under  the  conditions  of  a scale  height 
which  is  independent  of  altitude  and  winds  which  are  constant  in  time,  the  peak 
electron  density  steadily  decreases  to  its  asymptotic  value  and  never  increases. 
However,  in  a realistic  atmosphere  in  which  the  temperatures  and  the  mean 
molecular  weight  are  lower  at  lower  altitudes,  the  scale  height  decreases  at 
lower  altitudes.  The  peak  electron  density  can  then  be  expected  to  increase 
inversely  with  the  effective  number  density  scale  height. 

Altitude 'independent  but  time -varying  neutral  winds  do  not  have  a 
significant  effect  on  the  peak  electron  density . The  ion  cloud  just  evolves  to 
the  asymptotic  form  that  is  appropriate  for  that  altitude.  The  rate  at  which 
the  cloud  approaches  the  asymptotic  form  depends  on  the  altitude  at  which  the 
ion  cloud  finds  itself.  An  altitude -dependent  neutral  wind,  on  the  other  hand, 
caq^ause  changes  in  the  peak  electron  concentration.  An  altitude -dependent 
wind  at  lower  altitudes,  say  arcxind  150  km,  will  have  a greater  effect  on  the 
peak  electron  density  than  a wind  of  the  same  magnitude  at  a higher  altitude. 

The  change  in  the  peak  electron  density  would  be  gradual  and  not  sudden  for 
reasonable  models  of  the  neutral  wind.  We  should  note  that  even  in  Case  4 
for  a large  downward  wind  that  produced  the  highest  peak  electron  density  and 
narrowest  profile,  the  full -width  at  half-max  was  15  km  in  altitude. 

A nonvertical  magnetic  field  influences  the  behavior  of  small  ion  clouds 
more  than  large  ion  clouds.  However,  the  examples  we  have  shown  indicate 
that  for  large  dip  angles  and  vertical  neutral  winds,  there  is  little  difference 
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in  the  behavior  of  small,  low-conductivity  ion  clouds  compared  to 
high-conductivity  ion  clouds.  They  settle  under  the  influence  of  gravity  almost 
as  r^idly  and  achieve  essentially  identical  shapes  The  primary  difference 
between  low-density  and  high-density  ion  clouds  is  the  effect  that  electric  fields 
and  northward  neutral  winds  have  on  the  vertical  motion  of  the  ion  clouds.  As 
indicated  in  Eq.  (4.23)  both  westward  electric  fields  and  northward  neutral 
winds  cause  descent  of  the  ion  cloud  but  this  effect  is  less  pronounced  for  high- 
density  ion  clouds. 

As  the  cloud  changes  altitude,  the  field -line -integrated  Pedersen 
conductivity  of  the  cloud,  an  important  parameter  for  determining  its  motion 
and  striation  behavior,  changes.  For  our  constant  scale -height  atmospheric 
model  and  times  when  the  cloud  remains  at  sufficiently  high  altitude  so  that 
c = vM/eB  < < 1 , we  can  determine  the  time -dependence  of  the  conductivity 
explicitly.  By  neglecting  the  c in  the  denominator  in  Eq.  (4. 14),  making  use 
of  Eq.  (4.32)  for  u , and  denoting  the  initial  value  of  the  conductivity  when  the 
cloud  is  confined  to  an  altitude  range  small  compared  to  bv  2^  , we  may 
write  the  time -dependent  conductivity  as 
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In  deriving  the  above,  we  have  made  use  of  the  definitions  given  in  Eqs.  (4.36), 

1 ft 

(4.41),  (4.42),  (4.43),  (4.45),  (4.46),  (4.50)  and  (4.51)  and  the  evaluation 

/ 2\ 

2 f exp  |ir-j  dl  = t*  [1  + (v  + 1)  t*]  exp  . 

0 

This  remarkable  result  indicates  that  as  the  ion  cloud  is  evolving  and 
descending,  its  Pedersen  conductivity  is  given  solely  by  the  value  of  the 
altitude  of  the  center  of  the  ion  cloud.  If  there  is  a constant  downward  drift, 
the  Pedersen  conductivity  increases  exponentially . If  there  is  no  constant 
downward  drift,  the  Pedersen  conductivity  increases  linearly  with  time.  For 
the  values  of  the  parameters  chosen  here,  appropriate  for  the  Anne  cloud, 
the  doubling  time  for  the  conductivity  is  approximately  3/4  of  an  hour. 


5.  SPRUCE  ION  CLOUD  MODEL 


In  Section  2 we  provided  a description  of  the  e:q>ansion  of  the  neutral 
barium  atom  vapor  cloud.  In  Section  3 we  described  the  deposition  of 
ionization  on  magnetic  field  lines  resulting  from  the  photoionization  of  barium 
atoms  in  the  expanding  neutral  cloud.  In  Section  4 we  have  discussed  factors 
that  influenced  the  peak  electron  density  land  scale  size  of  the  ion  cloud 
parallel  to  the  magnetic  field  at  late  times. 

We  have  applied  the  procedures  developed  in  those  sections  to  obtain 
models  of  the  Spruce  ion  cloud.  We  will  first  list  the  parameters,  both  input 
and  derived,  associated  with  each  of  three  different  models.  A,  B and  C. 
Models  A and  B differ  only  in  their  late-time  properties.  Model  C,  on  the 
other  hand,  provides  what  we  believe  to  be  an  underestimate  of  the  total  ion 
inventory  in  the  Spruce  ion  cloud.  We  believe  that  the  characteristics  of  the 
Spruce  ion  cloud  are  best  exemplified  by  the  parameters  given  in  model  A 
(or  B).  Model  C is  presented  as  an  alternative  in  order  to  demonstrate  the 
sensitivity  of  the  models  to  some  of  the  input  assumptions. 

Table  5. 1 provides  a listing  of  the  parameters  associated  with  each  of 
the  three  models.  An  examination  of  the  first  ten  items  indicates  that  the 
input  parameters  at  early  time  for  models  A and  B are  identical.  Examination 
of  the  next  four  input  parameters  indicate  that  model  A and  B differ  only  in  the 
late-time  parameters.  Models  A and  C have  identical  late-time  input  para- 
meters. 

The  most  significant  input  parameter  is  item  5,  the  total  ion 
inventory.  Model  A results  from  the  attempt  to  maximize  the  total  ion 
inventory  and  still  have  consistency  with  other  reasonable  assumptions  and 
maintain  reasonable  agreement  between  the  derived  parameters  and  inter- 
pretations of  observations.  This  value  is  approximately  one -half  of  the 
inventory  assumed  in  our  previous  modeling.  ^ This  reduction  is  consistent 


*See  footnote  on  page  8. 
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Table  5. 1.  Parameters  for  Spruce  cloud  models 


Parameter 

(Units) 

Model 

A 

Model 

B 

Model 

C 

Inputs  to  Model  - early  time  {<  100s) 

1. 

Mass  of  barium  vapor,  Mo> 

(kg) 

3.56 

3.56 

2. 46 

2. 

Valorization  efficiency 

7.4% 

7. 4% 

5.1% 

3. 

Ionization  time  constant,  , 

(s) 

10 

10 

10 

4. 

Fraction  of  atoms  ionized,  f , 

(10^^ 

0.42 

0.42 

0.  42 

5. 

Total  ion  inventory 

6.  5 

6.  5 

4.5 

6. 

Atmospheric  density,  P- , 

(kg  Am®)  0.45 

0.45 

0.45 

7. 

Neutrjd  diffusion  coeffiment,  D,j , 

(km  Vs) 

0.052 

0.052 

0.052 

8. 

Equal  mass  radius,  a^ , 

(km) 

1.26 

1.26 

1.09 

9. 

Initial  expansion  velocity,  Rq  , 

(km/s) 

1.15 

1. 15 

1.20 

10. 

Onset  of  photoionization,  t. , 

(s) 

2.0 

2.0 

1.  5 

Inputs  to  Model  - late  time  (>100s) 

11. 

Parallel  diffusion  coefficient,  , 

(km^/s) 

0. 12 

0. 12 

0. 12 

12. 

Effective  scale  height,  , 

(km) 

30 

22 

30 

13. 

Initial  gravitational  drift,  gt^ , 

(m/s) 

8 

10.9 

8 

14. 

Downward  neutral  wind,  , 

(m/s) 

0 

8 

0 

Derived  Quantities 

1. 

Neutral  cloud  radius,  a,  at  t = t.  + T 

(km) 

3.  bo 

3.80 

3.51 

2. 

Neutral  cloud  radius  at  20  s 

(kn^ 

4.  72 

4.  72 

4.42 

3. 

Peak  ion  density  at  22  s 

(1^3  m- 

■3)  2.4 

2.4 

2.  2 

4. 

Maximum  ion  density 

(10^3  m' 

•3)  2. 45 

2.45 

2. 30 

5. 

Transverse  scale  size 

(km) 

2.  66 

2. 66 

2. 30 

6. 

Peak  ion  density  at  100s 

(10^3  HI- 

•3)  1.75 

1.75 

1.6 

7. 

Ion  cloud  length  at  100s 

Ml 

7.8 

7.8 

7.6 

8. 

Field  line  content 

(l^^m* 

■2)  2.4 

2.4 

2. 16 

9. 

Initial  Pedersen  conductivity 

(mh^ 

27.4 

27.4 

24.7 

10. 

Peak  ion  density  at  30  min. 

(10l3m' 

■3)  0.  G7 

0.70 

0.  51 

11. 

Ion  cloud  length  at  15  min. 

(km) 

18.5 

16.0 

18.5 

12. 

Ion  cloud  length  at  30  min. 

(km) 

24 

19 

24 
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with  our  present  modeling  and  a reduced  transverse  scale  size  for  the  ion 
cloud.  We  believe  that  the  modeling  presented  here  is  theoretically  more 
well-founded  and,  as  we  shall  see  below,  is  more  consistent  with  the  existing 
observations. 

However,  this  value  of  ion  inventory  is  considerably  larger  than  the 
24  21 

value  of  3.5  x 10  derived  by  Tsunoda  based  on  measurements  made  with 
the  incoherent -scatter  radar.  Since  we  will  be  making  a detailed  comparison 
of  the  results  of  models  A,  B and  C with  the  model  obtained  from  the 
incoherent-scatter  radar,  an  evaluation  of  the  reasons  for  the  different  values 
will  be  made  at  that  time.  The  ion  inventory  for  model  C was  arrived  at 
specifically  by  an  attempt  to  match  the  parameters  determined  by  the 
incoherent-scatter  radar.  We  believe  that  the  model  C value  represents  an 
actual  underestimate  of  the  total  ion  inventory  as  will  be  made  clearer  below. 

A unique  measurement  of  the  ion  inventory  in  the  Plum  ion  cloud 
produced  from  a different  48  kilogram  payload  barium  release  at  about  185 
kilometers  altitude  during  the  Secede  n test  series  was  made  by  Prenatt. 
Prenatt  made  dispersive -phase  measurements  by  receiving,  at  four  separate 
sites,  several  radio  frequencies  transmitted  by  a beacon  passing  behind  the 
ion  cloud.  These  measurements  provided  the  total  ion  content  along  four  cuts 
through  the  Plum  ion  cloud.  These  data  were  received  during  a time  span 
from  sq[)proximately  2 to  3 1/2  minutes  after  release.  The  Plum  ion  cloud  was 
already  quite  distorted  at  this  time,  but  Prenatt  was  able  to  successfully  model 
the  ionization  content  by  superposing  3 different  Gaussian  distributions. 

As  a result  of  his  modeling  Prenatt  determined  a characteristic  1/e 

dimension  transverse  to  the  magnetic  field  of  about  2.  8 km  and  a total  ion 

24 

inventory  of  7. 4 x 10  . Because  Plum  was  released  at  a slightly  lower 

altitude  than  Spruce  and  more  oxidation  of  Ba  by  the  larger  ambient  oxygen 

density  could  be  e3q)ected,  one  might  consider  this  value  of  ion  inventory  to  be 

a lower  limit  for  the  Spruce  ion  cloud.  It  was  the  existence  of  this  measure - 

24 

ment  that  motivated  us  to  adopt  a value  as  large  as  6. 5 x 10  for  the  ion 
inventory  in  model  A of  the  Spruce  ion  cloud. 
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The  value  of  the  ionization  time  constant  and  fraction  of  atoms 
ionized,  items  3 and  4 in  Table  5. 1,  were  discussed  in  Section  3.  We  have 
adopted  the  same  values  for  all  three  models.  Once  the  fraction  of  atoms 
ionized  is  assumed,  the  mass  of  the  released  barium  vapor,  M^,  is  deter> 
mined.  The  vaporization  efficiency,  item  2,  is  just  the  ratio  of  this  mass  to 
the  total  chemical  payload  weight.  ^ see  that  the  resulting  values  lie  within 
the  crude  5 to  10%  efficiency  estimates  that  are  often  discussed. 

The  atmospheric  density  model  adopted  determines  items  6 and  7. 

The  equa]-mass-radius,  a , is  defined  in  Section  2 and  derived  from  items  1 
and  6.  The  initial  expansion  velocity,  , affects  the  scale  size  of  the 
neutral  cloud  and  thus  the  ion  cloud.  Likewise  the  time  chosen  for  the  onset 
of  photoionization,  t^ , also  effects  the  scale  of  the  ion  cloud  and  the  sh^e 
of  the  density  distribution  transverse  to  the  magnetic  field.  In  Section  3 we 
indicated  that  a minimum  value  for  t^  was  1.  5 seconds.  A maximum  reason- 
able value  must  not  be  too  much  larger  than  2 seconds  because  radar 

observations  indicate  that  the  peak  electron  concentration  exceeds 
12  -3  9 

7 X 10  m within  3 seconds  of  release.  We  shall  defer  discussion  of  the 
late -time  inputs  to  the  models  until  we  discuss  parameters  derived  from 
these  inputs. 

We  now  discuss  the  results  obtained  by  applying  the  modeling 
concepts  developed  in  the  three  previous  sections  and  these  input  parameters 
to  the  Spruce  ion  cloud.  We  have  already  seen  in  Section  2 that  the  snowplow 
expansion  model  with  diffusion  added  provides  good  agreement  with  a photo- 
grj^h  of  the  neutral  cloud  at  20  seconds  after  release.  We  further  indicated 
that  diffusion  rapidly  causes  the  density  profile  to  become  Gaussian.  Further- 
more, the  cloud  radius  determined  by  the  snowplow  is  essentially  equal  to  a 
point  that  has  a density  1/e  of  the  central  density.  Hence,  we  approximate 
the  neutral  cloud  by  an  expanding  Gaussian  whose  Gaussian  radius  is  given  by 


a(t) 


0.63 


(5.1) 
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In  Section  3 we  described  the  deposition  of  ionization  on  magnetic 
field  lines.  An  appropriate  scaling  parameter  was  the  radius,  a , of  the 
neutral  cloud  at  t = tj.  + t.  , i.  e. , one  photoionization  time  after  photo- 
ionization is  assumed  to  begin.  The  resulting  values  of  a for  the  three 
models  are  listed  as  the  first  derived  quantity  in  Table  5. 1.  The  value  of 
this  radius  is  smaller  for  model  C because  the  mass  of  the  released  vapor 
is  smaller  and  the  time  at  which  photoionization  commences,  t.  , is  earlier. 
The  radius  of  the  neutral  cloud  at  20  seconds  after  release  is  listed  as  the 
second  derived  quantity. 

F^re  5. 1 shows  the  results  of  the  ion  deposition  calculation.  The 
four  solid  curves  show  the  total  ion  content  on  field  lines  as  a function  of  the 
position  across  the  magnetic  field  for  several  different  times  after  release. 
For  instance,  the  curve  labeled  t = 5 s represents  the  amount  of  ionization 
that  has  been  deposited  after  photoionization  has  been  operating  for  3 s.  The 
dashed  curve,  which  we  will  discuss  below,  is  a Gaussian  profile.  We 
observe,  as  remarked  earlier,  that  the  distribution  of  ionization  on  field 
lines  by  the  deposition  model  is  more  peaked  than  a Gaussian.  By  22 
seconds  after  release,  or  2 photoionization  times  t.  , the  total  content  is 
within  15%  of  the  total  content  deposited  on  field  lines. 

The  two  dashed  curves  in  Figure  5. 1 show  the  field-line  content  at 
two  separate  times  for  model  C.  The  total  ion  inventory  indicated  by  the 
dashed  curve  labeled  total  content  is  approximately  2/3  of  the  total  ion 
inventory  contained  within  the  solid  curve  labeled  total  content.  The  trans- 
verse scale  size  listed  as  the  fifth  derived  quantity  is  the  radius  at  which  the 
total  field  line  content  has  a value  of  1/e  of  the  field-line  content  passing 
through  the  center  of  the  cloud. 

Once  an  electron-ion  pair  is  created,  it  is  tied  to  the  magnetic  field 
line  on  which  it  originated.  It  can  no  longer  e}q}and  across  the  magnetic 
field  and  the  transverse  scale  size  of  the  ion  cloud  is  given  by  the  distribution 
shown  in  Figure  5. 1.  There  is  no  such  restriction  on  the  motion  of  the  ions 
parallel  to  the  magnetic  field.  During  the  first  20  seconds  after  release  we 


(0 

H < O 

|§  8 

8 s s 


a ww j j j 
g g^e  a a 

si^ii  i 


r 


/I',  i 


(jW/SNOi  ^jOi)  iNaxNOO  ann  oiaia 


SNOI  ,.0I)  • ((  ^ XV  AXlSNSa  NOI 


V U O’ 

'm  i "H  :§  ® 

g“”5”S 

I S'®  s is 

S n §!  ®|  «- 
3 oj  8^ 

n o u ® ^ m 

gfi-sIsS 

3 t<  X3  ® g ^ 

TJ  ® _ "O 

S m -ot®  S c 
o « c ^ o 
W ® a o p c 

js  t-  iS  ® ^ j3 

it  S fH  T3  n 
. ® ® ® ® o 
O *3  ® ^ ®‘ 

^ mH  §-§5 
a ® t-i  t2 

CQ'Sr®  si  I 

r c o 2 ® ^ 

<5  ® ■*■*  o •>  M 

oj  b ® >H  ® ® 

■o  S ® ^ o _ 
0^3  > A-  W '3 

B'S'i  e ®*® 

o ^ P -o 

JS  ■&  lU  *0  33  M 

K«'~'*'  - ® ® . 

^ ^ g5|^  § 

■s-s  11 

a *2  ^ ® x;  ®18 
® g o bD  ^ Eh  ^ 
g g g 


g 'S  § S -O  .:§ 

Jig’S 

rt  ® O ®’^ 

C T3  --  o tj  m 

•S-S  a 2 ® I ® 

•a  ^ cB  >* 

w flj  G ® f?  2 ^ 

•o  ^ "g  M g?  Stj 

g jg  a d g 1 2 

»2  m O-  o H u ts 


assume  that  the  ion  distribution  parallel  to  the  magnetic  field  is  the  same  as 
that  given  by  the  neutral  cloud  from  which  the  ionization  originated.  Knowing 
the  content  on  the  field  line  and  the  scale  length  parallel  to  the  field  given  by 
Eq.  (5.1)  we  can  derive  an  ion  concentration  profile  at  approximately  22  sec- 
onds after  release.  The  left-hand  scale  for  the  ion  concentration  should  be 
read  only  for  the  solid  and  dashed  curves  labeled  t = 22  s. 

Also  shown  in  Figure  5 . 1 by  a dotted  curve  is  the  Gaussian  model 
21 

derived  by  Tsunoda  based  on  measurements  with  the  incoherent-scatter 
radar.  This  dotted  curve  can  be  compared  to  the  solid  and  dashed  curves 
labeled  t = 22  s.  We  see  that  the  model  C peak  ion  concentration  agrees 
with  the  Tsunoda  model  value,  but  that  its  Gaussian  dimension  trans- 
verse to  B is  a bit  small.  Model  A,  on  the  other  hand,  has  a peak  ion  con- 
centration that  exceeds  the  Tsunoda  model  value  by  ^proximately  10%  but  pro- 
vides an  excellent  agreement  with  the  measured  scale  size  transverse  to  S . 

If  the  ion  cloud  actually  had  the  shape  given  by  model  C,  the  incoherent-scatter 
radar  would  probably  have  measured  a lower  peak  ion  concentration  due  to  the 
finite  beam  width.  It  is  for  this  reason  that  we  believe  model  C represents  an 

imderestimate  to  the  actual  size  of  the  ion  cloud  and  that  the  actual  ion  inven- 
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tory  is  in  excess  of  the  model  C ion  inventory  of  4.5  x 10  ions. 

Ionization  is  85%  complete  by  22  seconds  and  we  assume  that  the  ion 
cloud  expands  parallel  to  the  magnetic  field  by  diffusion,  i.e. , we  assume  that 
the  expansion  due  to  the  snowplow  ceases.  The  increase  in  length  of  the  ion 
cloud  along  the  magnetic  field  is  shown  by  the  dotted  line  in  Figure  2.4.  We 
note  that  this  increase  is  not  substantially  different  than  that  predicted  by  the 
snowplow  without  diffusion.  Thus  initially  we  can  represent  the  increase  in 
the  Gaussian  scale  size  of  the  ion  cloud  parallel  to  the  magnetic  field  by 


The  time -dependence  of  the  total  ion  inventory  is 


N{t)  = N^[l  -exp(-(t-tj)/T.)] 


(5.3) 


for  t > t.  . We  may  develop  a model  for  the  time  dependence  of  the  peak  ion 
concentration,  n^  , initially  by  taking  into  account  the  increase  due  to  the 
additional,  remaining  ionization,  and  the  decrease  due  to  diffusion  parallel  to 
the  magnetic  field.  We  obtain  an  expression  of  the  form 

from  which  the  maximum  peak  ion  density,  listed  as  item  4 in  Table  5.1,  has 
been  determined. 

We  have  examined  the  evolution  of  the  Spruce  ion  cloud  as  it  settles 
under  gravity.  Figure  5.2  shows  the  development  of  the  Spruce  ion  cloud 
according  to  model  A.  Initially,  for  the  first  10  minutes,  the  peak  ion  con- 
centration behaves  according  to  Eq.  (5.4).  This  equation  is  a valid  description 
of  the  time  dependence  of  the  peak  concentration  because  the  ion  cloud  has  not 
descended  very  far  into  the  atmosphere  and  the  diffusion  coefficient,  , has 
essentially  its  initial  value  given  by  the  input  parameter  item  11  in  Table  5.1. 
The  peak  ion  concentration  and  the  ion  cloud  length  at  100  seconds,  listed  as 
items  6 and  7 in  Table  5.1,  are  given  by  Eqs.  (5.4)  and  (5.  2)  respectively.  At 
later  time,  however,  when  the  parallel  diffusion  coefficient  evaluated  at  the 
center  of  the  ion  cloud  is  different  from  its  initial  value,  the  peak  ion  concen- 
tration ceases  to  fall  as  rapidly  as  predicted  by  Eq.  (5.4). 

Figure  5 . 3 shows  a plot  of  the  peak  ion  concentration  as  a function  of 

time  after  release.  The  late -time  curve  for  model  A results  from  the 

ion  cloud  settling  imder  the  influence  of  gravity.  We  note  that  the  ion  concen- 
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tration  has  fallen  below  7 x 10  m by  18  minutes  alter  release  for  model  A. 

90 

The  short-pulse  radar  operated  by  Oetzel  and  Chang  continued  receiving 
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1 0»  • 1 3 
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ION  density  (n«»-3) 

Figure  5. 2.  Expansion  and  descent  of  the  Spruce  ion  cloud  for  model  A under 
the  assumption  of  no  winds  and  constant  effective  neutral  density 
scale  height.  Times  shown  are  times  after  release. 
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returns,  although  sporadically,  at  24  MHz  indicating  the  presence  of  an  over- 
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dense  cloud  with  an  ion  concentration  of  7 x 10  m as  late  as  30  minutes 
after  release. 

Since  downward  winds,  or  a descent  of  the  cloud  to  a region  of  cooler 

atmosphere  where  the  atmospheric  scale-height  is  less,  decrease  the  effective 

diffusion  coefficient,  the  ion  cloud  would  not  expand  along  the  field  line  as  rapidly 

as  predicted  by  model  A.  In  order  to  model  this  effect  we  adjusted  the  late-time 

input  parameters  for  model  B.  We  decreased  the  neutral  density  scale -height 

to  22  km  (the  same  value  used  in  Section  4 for  the  Anne  ion  cloud)  and  added  a 

downward  neutral  wind  of  8 m/s.  In  order  to  keep  the  initial  diffusion  coefficient 

the  same  in  this  cooler  atmospheric  model,  we  had  to  reduce  the  atmospheric 

density  at  190  km  by  a small  amount  resulting  in  the  increased  value  of  g • 

By  making  these  changes  we  find  that  the  model  B curve  does  maintain  the  ion 
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concentration  above  7 x 10  m for  more  than  30  minutes. 

We  note  that  the  model  A and  B peak  ion  concentration  curves  pass 

above  the  measurements  made  at  22  seconds  and  100  seconds  after  release  by 

24 

two  Independent  radar  systems.  The  HF  long -pulse -radar  value  does  provide 

only  a lower  bound  to  the  peak  concentration.  Other  48  kg  barium  releases 

have  indicated  peak  values  at  100  seconds  after  release  of  between  1 . 6 and 
13  -3 

1.9x10  m .It  was  primarily  for  this  reason  that  the  ion  inventory  for 
model  A was  chosen  so  that  the  model  A curve  would  lie  slightly  above  the 
measurement  made  at  100  s. 

In  model  C the  inventory  was  reduced  so  that  the  peak  ion  concentration 
would  agree  with  these  two  early -time  measurements.  We  have  already  indi- 
cated when  we  discussed  Figure  5 . 1 why  we  feel  that  the  density  profile  predicted 
by  model  C would  actually  give  a measured  value  at  22  seconds  less  than  the 
value  that  was  measured.  The  model  C cuive  essentially  lies  10%  below  the 
model  A curve. 

Model  B not  only  predicts  a less  r^id  decrease  in  the  peak  ion 
concentration  but  it  also  predicts  that  the  ion  cloud  Gaussian  length  parallel  to 


i 
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the  magnetic  field  does  not  increase  so  rapidly . The  two  solid  curves  in 
Figure  5.4  show  the  increase  in  Gaussian  length  as  a function  of  time  for 
models  A and  B respectively.  By  comparison,  we  have  shown  the  increase 
in  length  predicted  by  Eq.  (5 . 2)  for  two  different  values  of  the  parallel  diffusion 
coefficient  as  indicated  by  the  two  dashed  curves.  We  note  that  for  the  first  ten 
minutes  after  release  there  is  little  difference  between  the  four  curves.  How- 
ever, measurements  of  the  full  width  of  the  ion  cloud  parallel  to  the  magnetic 
field  using  the  incoherent-scatter  radar  have  indicated  that  by  about  15  to  20 
minutes  after  release  the  Gaussian  length  was  16  to  18  km  and  did  not  increase 
very  rapidly  after  this  time . This  measurement,  in  addition  to  the  short -pulse 
radar  measurement,  lends  support  to  the  model  B description  of  the  late -time 
behavior  of  the  Spruce  ion  cloud. 

The  field-line -content  for  the  three  models  is  given  as  item  8 in  Table 
5.1.  The  initial  field -line -integrated  Pedersen  conductivity,  defined  by  Eq.  (4.14), 
can  be  calculated  from  a knowledge  of  the  field -line -content  and  is  given  as  item 
9 of  Table  5.1.  At  later  times,  as  the  ion  cloud  descends,  the  Pedersen  conduc- 
tivity increases  according  to  Eq.  (4.57).  The  last  three  items  in  Table  5.1  list 
some  of  the  late -time  parameters  derived  from  these  three  models. 

The  procedures  that  have  been  developed  in  Sections  2,  3 and  4 have 
provided  a satisfactory  explanation  of  many  of  the  observations  of  the  Spruce 
ion  cloud.  These  procedures  can  be  extended  to  model  other  ion  clouds.  The 
real  test  of  the  modeling  concepts  that  have  been  developed  here  is  to  apply  the 
model  to  ion  clouds  of  different  payloads  released  at  different  altitudes.  So  far 
this  task  has  not  been  done . However,  based  on  the  success  that  we  have 
achieved  with  this  model,  we  are  confident  in  applying  it  to  other  ion  clouds  of 
around  48  kg  chemical  payload  weight  released  in  the  vicinity  of  180  to  190  km 
altitude. 
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Figure  5. 4.  Spruce  ion  cloud  Gaussian  length  parallel  to  the  magnetic  field 
as  a function  of  time  after  release  for  models  A and  B as 
compared  to  the  assumption  of  constant  parallel  diffusion 
coefficients. 
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